Geometry-dependent critical currents in superconducting nanocircuits 
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In this paper we calculate the critical currents in thin superconducting strips with sharp right- 
angle turns, 180-dcgrcc turnarounds, and more complicated geometries, where all the line widths 
are much smaller than the Pearl length A = 2A^/d. We define the critical current as the current 
that reduces the Gibbs free-energy barrier to zero. We show that current crowding, which occurs 
whenever the current rounds a sharp turn, tends to reduce the critical current, but we also show that 
when the radius of curvature is less than the coherence length this effect is partially compensated 
by a radius-of-curvature effect. We propose several patterns with rounded corners to avoid critical- 
current reduction due to current crowding. These results are relevant to superconducting nanowire 
single-photon detectors, where they suggest a means of improving the bias conditions and reducing 
dark counts. These results also have relevance to normal-metal nanocircuits, as these patterns can 
reduce the electrical resistance, electromigration, and hot spots caused by nonuniform heating. 
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I. INTRODUCTION 

When an electrical current travels through a 180- 
degree hairpin turn or around a sharp corner in a thin 
film, it tends to concentrate on the inner boundary of 
the curve. For normal metals, this effect is known as 
current crowding, and results in an additional resistance 
over what would be expected from naive geometrical 
arguments.^ The situation with a superconductor is anal- 
ogous to that in normal metals, except that the effect 
manifests itself primarily as a reduced critical current for 
the pattern. A superconducting wire will become resis- 
tive as soon as the critical current for the inner corner 
is exceeded, even if the total current is lower than the 
critical current of the connecting straight-line segments. 

Current crowding has an important implication for 
measurements of critical currents in thin films because 
such measurements typically include at least one sharp 
corner. For example, in two-point measurements, nar- 
row, straight segments connect to larger contacts at the 
end, forming sharp inner corners on either side. Even the 
ubiquitous Kelvin Bridge resistor striicture typically in- 
cludes sharp corners where the voltage leads contact the 
center current-carrying lead. These measurements are 
used extensively in materials and device characterization 
for nanoscale structures, but in practice these measure- 
ment methods can in general underestimate the critical 
current. 

Nearly 50 years ago, Hagedorn and Hall^ extended 
earlier work by Cockroft^ to consider a strip conductor 
with a right-angle bend. They calculated both the re- 
sistance for a normal-metal strip and the current-density 

distribution for a superconducting strip. Using confor- 
mal transformations, the authors showed that the cur- 



rent density increases on the inside corner of a bend, and 
derived the mathematical form for an optimally rounded 
inner boundary for which the critical current for the bend 
would be the same as for a straight strip. The authors' 
resistance result and its extension to other geometries^ 
became widely used in the field of integrated-circuit lay- 
out, but their critical-current calculations were appar- 
ently forgotten. Despite extensive work on the critical 
current of superconducting materials over the interven- 
ing period,^ to our knowledge Hagedorn and Hall's pre- 
dictions have never been tested experimentally. 

Although Hagedorn and Hall correctly predicted the 
presence of current crowding, their treatment did not 
permit a quantitative estimate of its effect on the crit- 
ical current, and was limited to considering only simple 
90° corner goemctries. A new aspect emphasized in our 
paper, not considered in Rcf. 1, is that the critical cur- 
rent of a superconducting strip in the absence of thermal 
excitations is the current at which a nucleating vortex 
surmounts the Gibbs free-energy barrier at the wire edge 
and then is driven entirely across the strip. The result- 
ing voltage along the strip is proportional to the rate at 
which these processes occur. Our approach permits an 
estimate of the critical current of a variety of thin-film 
patterns. To calculate the critical current, one needs to 
know both the current distribution and the vortex's inter- 
action energy with the edge of the strip. An interesting 
and important consequence of this treatment is that the 
critical current of a right-angle bend is finite, even though 
the current density at the sharp inner corner calculated 
using the London equation, as done in the Hagedorn pa- 
per, would diverge. In this paper, we perform theoret- 
ical calculations of the critical current using conformal 
transformations going beyond Hagedorn and Hall's cal- 
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culations for simple right-angle bends to include many 
more pattern geometries. 

Thermal excitation of vortices over the Gibbs free- 
energy barrier can produce random voltage pulses due 
to vortex motion across the strip. ^ The resulting time- 
averaged voltage gives rise to a strongly temperature- 
dependent finite resistance of the superconductor at cur- 
rents below the critical current that would occur in 
the absc;nc;c; of thermal excitations. Voltage pulses like 
these have been proposed as the origin of so-called dark 
counts in superconducting nanowire single-photon de- 
tectors (SNSPDs).^"^^ Our results should provide guid- 
ance regarding the extent to which thermally activated 
voltages depend upon the geometry of the underlying 
nanocircuitry. 

These results may have important implications for 
SNSPD performance. SNSPDs consist of thin films pat- 
terned to form a continuous wire arranged in a boustro- 
phedonic pattern on the substrate. Their performance 
depends sensitively on the bias current: Typically, de- 
tectors increase in sensitivity with increased bias current 
until the critical current is exceeded. Commonly, some- 
where between a few and 100% of fabricated devices suf- 
fer from an effect known as "constriction," in which the 
critical current is anomalously suppressed. This suppres- 
sion is established by comparing the measured critical 
current to one predicted by measuring the device induc- 
tance and fitting to theoretical expectations (which can 
be used to extract a "true" critical current). At present, 
the cause of constrictions in SNSPDs is not understood. 
There are no unambiguous demonstrations in which a 
constriction was identified. It seems from our work that 
sharp corners, hairpin turns, contact inner corners, and 
line asperities (both widenings and narrowings) could all 
result in constriction. Our results, as well as the earlier 
results of Hagedorn and Hall, are of direct relevance not 
only to superconducting wires, but also to conventional 
current flow in patterned thin films. Sharp corners and 
current crowding are known to cause problems of electro- 
migration and the formation of hot spots in the pattern. 
Such effects can limit the performance and reliability of 
integrated circuits. The results shown here specify how 
to optimally design a variety of patterns without current 
crowding, thus minimizing the undesirable impact of this 
effect. 

In Sec. II, we show why the current density is very 
nearly uniform in straight, narrow, thin superconducting 
strips. In Sec. Ill, we define the critical current as the 
current that lowers to zero the Gibbs free energy bar- 
rier against vortex entry into the strip. In Sec. IV, we 
present a general procedure for calculating the critical 
current and, as an example, apply it to the case of a long, 
straight, narrow, thin superconducting strip. In Sec. V, 
we show how current crowding, which occurs whenever 
the current bends around a curve, tends to reduce the 
critical current. In Sec. VI, we show how the current- 
crowding effect is ameliorated by a radius-of-curvature 
effect. In Sec. VII, we calculate the extent to which the 



critical current is reduced in 180-degree turnaroimds. We 
show that sharp inner corners reduce the critical cur- 
rent, and we show how the inner curve can be optimally 
rounded so that there is no critical-current reduction. 
In Sec. VIII, we calculate the extent to which the criti- 
cal current is reduced in 90-degree turns. We show that 
sharp inner corners reduce the critical current, and we 
show how the inner curve can be optimally rounded, as 
noted in Ref. 1, so that there is no critical-current re- 
duction. In Sec. IX, we show how current crowding at 
the sharp corners of sidebar contacts in the shape of T 
intersections reduce the critical current, and in Sec. X, 
we show how a similar critical-current reduction occurs 
because of ciirrent crowding at the sharp corner where 
a long strip makes electrical contact with a wide con- 
tact pad. In both Sees. IX and X, we also suggest new 
patterns with rounded corners that should avoid critical- 
current reduction. In Sec. XI, we examine the role played 
by edge defects in reducing the critical current. In Sec. 
XII, we show how thermal excitation over the Gibbs 
free-energy barrier leads to broadening of the switching- 
current distribution observed when the critical current is 
measured repeatedly. In Sec. XIII, we compare our calcu- 
lations with prcvioiisly unexplained experimental results 
by Yang et al.^"' In Sec. XIV, we list some limitations 
of our results and discuss how these might be improved, 
and we examine some practical consequences. Appendix 
A contains a calculation of the self-energy of a Pearl vor- 
tex in a sector of arbitrary angle a. 



II. NEARLY UNIFORM CURRENT DENSITY 
IN A STRAIGHT, NARROW, THIN 
SUPERCONDUCTING STRIP 

We focus our attention on the properties of a thin su- 
perconducting strip of thickness d much less than the 
London penetration depth A, width W much less than 
the Pearl length^^ A = /d, and Ginzburg-Landau co- 
herence length ^ miich less than the width. In particular, 
let us examine the current density in an infinitely long 
thin superconducting strip of uniform width W in the 
xy plane, centered on the x axis. Suppose a current / 
flows in the x direction. Since the film thickness obeys 
d <^ \, the current density j is very nearly independent 
of z across the film thickness, and therefore it is conve- 
nient to consider only the sheet current (the integral of j 
across the thickness) K = jd. When W ^ K, K \s very 
nearly independent of y {\y\ < W/2). 

According to the London equation, the sheet- 
current density in the strip obeys K = — (2//ioA)[A -|- 
(0o/27r)V7], where A is the vector potential {B = 
V X ^4) and 7 is the phase of the order parameter. 
When the strip contains no vortices but carries a cur- 
rent / = /l^^^/2 ^x{y)dy in the x direction, the gauge of 
the vector potential can be chosen to absorb the constant 
((!i>o/27r)V7, such that Kr^(y) = — 2Aj.(y)//ioA, which is 
an even function and the self-field Bz{y) = —dA^/dy is 
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an odd function of y. That Kx is very nearly independent 
of y when ^ A can be shown using perturbation the- 
ory by expanding A',. = -2Ax/noh = Kxq + Kxi+Kx2 + 
... and Ax = A^o + ^xi + + where K^o = I /W , 
Kxn = — 2>la;„(y)//ioA is proportional to (W^/A)", and 

!^w/2^^n{y)dy = -2j^^^^^Axn{y)dy/no^ = for 
n > 1. To obtain the first-order correction term Kxi, 
we use the Biot-Savart law with Kxo to obtain 



HoKxo , 
'B^i(y) = ^^In 



y + W/2 



y-w/2 



(1) 



and then integrate Szi(y) = —dAxi/dy with respect to y 

for \y\ < W/2 to obtain Axiiy). The resuh is = 
-2Axi{y)/iio^ = {W/A)f{2y/W)Kxo, where 

f{u) = ;^[l-ln4+(l+u)ln(l+u)+(l-M)ln(l-u)]. (2) 

The function f{2y/W) varies smoothly from a minimum 
of /(O) = -0.061 at the middle of the strip to /(±1) = 
0.159 at the edges of the strip. Thus for the experiments 
reported in this paper, for which W < 100 nm and A « 
100 /um, the current density Kx{y) vs y in a long, straight 
superconducting strip is flat to better than 0.02%. 

In general, the current flow in the strip obeys V-JC = 
and (via the London equation) V x If = —2B/noA. 
However, as we have seen above, when W <^ A, the 
effect of B on the current flow is negligible, and to ex- 
cellent approximation the current flow can be obtained 
from V ■ K = Q and V x if = 0, as in thin non- 
superconducting strips. As has been demonstrated in 
Refs. 1-3, the current distributions in nonsuperconduct- 
ing thin films with complicated geometries can be ob- 
tained using complex-field methods. Here we express 
the X and y components of the sheet-current density 
K{x,y) = xKx{x,y) + yKy{x,y) as the real and imag- 
inary parts of a complex sheet-current density 1C{Q = 
Kx{x,y) — iKy{x,y), where /C(C) is an analytic function 
of the complex variable C = x+iy. The Cauchy-Riemann 
conditions obeyed by analytic functions correspond to the 
physical conditions that V ■ K = and V x K = 0. 



III. PROCEDURE FOR ESTIMATING THE 
CRITICAL CURRENT OF NARROW 
SUPERCONDUCTING STRIPS 

Within the context of the Ginzburg-Landau (GL) 
theory,^^"^° the upper limit to the dc supercurrent 
that can be carried before a voltage appears along the 
length of the sample is the GL critical current den- 
sity, Jgl = 'J^'o/S^^^TT^oC^^, critical sheet-current den- 
sity, Kgl = 2(/)o/3^/^7r/ioC-^j or critical current, Igl = 
2(/)oVr/3'^/^7r/io^A. However, since Jql was calculated 
theoretically assuming that both d and W were much 
smaller than both A and ^, these results cannot be ap- 
plied to the case of interest here, ^ <^ W. 

In this paper we use a general but approximate method 
for estimating the critical current at which a voltage 



first appears along the sample length, the voltage be- 
ing produced by the micleation of a vortex at an edge 
or corner and its siibsequent transit across the width. 
The method must be capable of calculating the criti- 
cal current for samples of many different geometries, in- 
cluding width variations and turnarounds. We will start 
by writing down the Gibbs free energy G(r^), where r„ 
is the two-dimensional coordinate of a vortex, G{ry) = 
Eseijiry) - Wi{ry). Eseif{ry) is the self-energy of the 
vortex accounting for all its interactions with the sample 
edges, including image vortices, and Wi{ry) is the work 
done by the sources of the current in moving the vortex 
away from the sample edge to the position r^. 

Neglecting core contributions, we use the London 
model of a vortex to calculate the self-energy Egei / as the 
area integral of the kinetic-energy density per unit area 
/ioAif^/4 outside the vortex core, taken to have radius 
^, the Ginzburg-Landau coherence length. ^'^'^^^^^"^'^ The 
divergence theorem can be used to obtain Es,,if{ry) = 
4>alcirc{i'v)l'^i where Icirc{i"v) is the net self-generated su- 
percurrent circulating around the vortex core when the 
vortex is at r„. 

When the vortex is placed in a sample carrying a cur- 
rent / distributed as the sheet-current density Kj{r), 
the work W7(r„) is the line integral of the Lorentz force 
FLif) = Kj(r) X (f)()Z from the vortex's point of entry to 
the position r„. Thus VK/(r„) = (pQAI{ry), where A/(r„) 
is the portion of the current / that fiows between the vor- 
tex's point of entry and the position r„. In contrast to 
the approach used in Refs. 10 and 11, we neglect any 
suppression of the magnitude of the GL order parameter 
by the current density. 

We wish to be able to calculate the critical current for 
a variety of sample geometries. In each case wc will find 
that for small currents / there is a Gibbs free-energy bar- 
rier at rb that prevents the nucleation of a vortex. At the 
top of the barrier, where = ri,, VG{ry) = 0, which cor- 
responds to a balance of forces: The Lorentz force, which 
tends to repel the vortex away from the sample edge, is 
balanced by the image force, which tends to attract the 
vortex back to the edge. We will define the critical cur- 
rent Ic as that value of the current for which the height 
of the barrier is reduced to zero and r^, moves to the po- 
sition Tc, where G{rc) = 0. In the following calculations 
we will find for each case that r,. is only a short distance, 
of the order of ^, from the vortex entry point. Since we 
consider ^ <^ W, the functions needed to calculate G{ry) 
can all be obtained using power-law expansions. 



IV. CRITICAL CURRENT DENSITY OF A 
LONG, STRAIGHT, NARROW, THIN 
SUPERCONDUCTING STRIP 

In this section, we present a general procedure for cal- 
culating the critical current in all geometries but apply 

it as an example to the simplest case, a long, thin strip 
of uniform width W , here assumed to occupy the space 
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< y < VF.^^ In general, we will first use conformal map- 
ping to find the applied sheet-current distribution and to 
calculate Wi{ry) for vortex positions close to the nu- 
cleation point. Next we will use the same mapping to 
find the current distribution around the nucleating vor- 
tex and to calculate Egeifivv). We will then examine the 
Gibbs free energy G(r„) for r„ at or near the barrier, 
and we will define the critical current as the current that 
reduces the Gibbs frcc-cncrgy barrier to zero. 
The conformal mapping, 



Ciw) 

w{C) 



dw 



2W 



<w + 



7r(u>2 — 1) ' 



(3) 
(4) 
(5) 



maps points in the upper half w-plane v > {w = u + iv) 
into the strip <y <W in the ^-plane = x + iy). 
The complex potential 



I fw — 1 
Qwiw) = - In — 

TT \ W + 1 



(6) 



describes the flow of current I in the upper half w-plane 
from a source at t« = 1 to a drain at w = — 1. In the 
C-plane the same complex potential is 



GdC) = Gw{w{0) 



TT Vu;(C)-|-l. 



(7) 



and the corresponding complex sheet current is /C(^(C) = 
dGc{C)/d( = Kx{x,y) — iKy{x,y). The streamlines of 
the applied sheet-current density K = xK^ + yKy arc 
obtained as contours of the stream function S{x, y) = 
^Qi^{x + iy), the imaginary part of 5^(C)- In general, 
KJix,y) = dS{x,y)/dy and Ky{x,y) = -dS(x,y)/dx. 

In the example for which w{C,) is given by Eq. (5), 
we have, with / = KiW, 5^(0 = KiC, IC^^lc) = Ki, 
S{x, y) = Kjy, = Kj, and Ky = 0. 

The complex potential describing the sheet-current 
flowing aroimd a Pearl vortex"'^^ sX w = in the up- 
per half ?y;-plane, subject to the boundary condition that 
at i; = there be no current flow perpendicular to the 
u axis, can be obtained simply by the method of im- 
ages, so long as the vortex's distance from the edge of 
the fllm is much less than the Pearl length A.^^ The self- 
generated sheet-current density circulating at distances 
p A aroimd the axis of a Pearl vortex has magnitude 
K = (j)o/nfj,oAp}^ We therefore have 



Gvw{wv;w) = 



»0O 
TT^oA 



In 



(8) 



where the term in the numerator arises from the nega- 
tive image at w = w* needed to satisfy the boundary 
condition at w = 0. The complex potential describ- 
ing the sheet-current flow circulating around a vortex at 



Cv = Xv + iyv in the C plane can be obtained by starting 
from Eq. (8) and using w{() from Eq. (5) to obtain 



icpo 



In 



TT/ioA V w{C) - w{Cy) 



)■ 



(9) 



This function automatically accounts for the infinite sets 

of positive and negative image vortices needed to satisfy 
the boundary condition that the sheet-current density be 
parallel to the edges at y = and y = W. The imagi- 
nary part of Qy(^ is the stream function Sy{xy, y^: x, y) = 
'^GvciCv'X)^ which rises to its largest value as ^ — > C?; 
and w{() — >■ w(Cj;); it is zero on the boundaries y = 
and y = W. 

The current circulating around the vortex can be ob- 
tained from Eq. (9) by evaluating the stream function 
at a cutoff radius equal to the Ginzburg-Landau coher- 
ence length, i.e., a distance ^ from the vortex position 

Ct,, where we make the approximation that ^ <C |Cu|- In 
the numerator of the argument of the logarithm we re- 
place w{Q by w{(y), and in the denominator we replace 

w{0 - wiC,) by ^du;(C„)/dC. = e/C'(«^(C))- When the 
vortex is at {x,y) — {x^,y.u), the circulating current is 



"^circip^ViUv^ — SviXv^yv-} Xv ^ ^lyv) 



-In 



k(C.)-«^*(C.)IIC'(^(C.))|- 



.(10) 



We can use the stream function S{x,y) to calculate 
the work Wi{ry) = (j)oAI{ry) done by the source of 
the current in moving the vortex from its point of en- 
try {x,y) = {xen,yen) to the point {x,y) = {xy,yy), 



AI{Xy, yy) — S{Xen: Ven) ~ S{Xy,yy). 



(11) 



Note that A/ also can be obtained by integrating 
the sheet-current density Ki{r) that passes between 
{x en, Ven) and (a;„,?/„). 

When ^ <C W . the distance between the point of entry 
and the point where the Gibbs free energy is zero is only 
a little larger than ^ for all the cases considered here. 
We therefore will use expansions of w{C) valid at small 
distances from the vortex entry point to evaluate Eqs. 
(7)-(ll) for all cases considered in this paper. 

As an example, let us now apply the above general 
procedure to evaluate the critical current for a long 
straight strip of width W . If the vortex entry point is 
at {xen, Ven) = (0, W), which Corresponds to C, = iW and 
w = 0, and the vortex position of interest is (xt,,y„) = 
(0, W — 5), wc can use Eq. (3) to obtain the approxima- 
tion w = m5/2W, which is valid for 5 <C W. Using this 
w in place of w{C,) in Eqs. (7) and (9) and evaluating Eqs. 
(10) and (11) by expanding them to lowest order in 5/W 
yields 



G 



In 



'25\ 



(12) 



The first term on the right-hand side of Eq. (12) is the 
self-energy accounting for the vortex's interaction with 
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its nearest negative image vortex at a distance 2S, and 
the second term is the negative of the work done by the 
source of the current in moving the vortex a distance S 
in from the edge. 

The free energy barrier occurs at (5 = 6b. Setting 
dG/d6 = there, we obtain 



(13) 



which is equivalent to the force balance between the re- 
pulsive Lorentz force ^o-?^/ and the attractive force of 
the image vortex (/)o/27r/LioAs. Setting G = at yields 
5i) = 5c = e^/2 = 1.36^ and the critical sheet current 
Ki = Kc, 




(14) 



where e is Euler's number, 2.718.... Despite the ap- 
proximations made in deriving Eq. (14), this result for 
Kc is numerically close to Kgl- The barrier height for 

Ki < Ka is 



In 



(15) 



In later sections we have followed the above proce- 
dure, employing complex fields and conformal mapping, 
to calculate the current flow and the critical sheet current 
in strips with various geometries, including turns and 
turnarounds. The main difference between these cascis is 
the mathematical form of the conformal mappings, which 
replace Eqs. (3)- (5). 
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FIG. 1: Current flow in a strip carrying average sheet-current 
density Ki around a circular arc of inner radius a and outer 
radius h = 2a, showing current crowding near the inner radius. 



a nucleating vortex from the point {x,y) = (a,0), the 
self-energy of the vortex, including its interaction with 
the nearest negative image, is the same as the first term 
in Eq. (12) except for correction terms of order d/a, which 
are negligible in critical-current calculations when ^ <C a. 
Thus when S <^ {b — a), the Gibbs free energy becomes 



G 



27r/xoA V ^ 



(y) - 0oi^0(a)(5. 



(16) 



Following the steps that led to Eq. (14), we obtain 



V. CURRENT CROWDING 

When a superconducting strip of constant width W 
follows a curving path, current crowding occurs near the 
inner radius of a bend. As a result, the critical current of 
the strip is reduced because vortices more readily nucle- 
ate where the edge has the minimum radius of curvature. 
Here we examine this efii'ect in two limits, ^ <^ W and 



A. London-model calculation for ^ <^W 

This critical-current reduction can be understood us- 
ing the London model for the geometry shown in Fig. 1. 
When a shcet-ciuTcnt K(p) = K^{p)(f> flows around an 
annulus of inner and outer radii a and b = a+W, where p 
is the radial coordinate and W <^ A, the solution obeying 
V ■ K ^ and V X K = is K^{p) = {b~a)Ki/ \n{b/a)p, 
where Kj is the average sheet-current density. The cur- 
rent density at the inner radius, whose radius of cur- 
vature is pc = a, has the enhanced value, K^[a) = 
{b — a)Kj/lii{b/a)a. As a function of the distance 6 of 



and 5c = eC/2 = 1.36^, (17) 



which tells us that when <^ pc = a, the critical current 
is reached when K^{a) = (^o/e7r/zoCA. However, since 
Ki < K^{a) and Kc is the critical value of Ki, we have 



Kr = 



eTT^oCA 



R, where R = 



aln(6/a) 
{h-a) 



(18) 



is the reduction factor due to current crowding at the 
curving inner boundary, relative to the Kc for a long 
straight strip [Eq. (14)]. For example, R = 0.69 when 
(& — a) = W = a. Note that i? — >■ 1 in the limit of infinite 
radius of curvature {pc = a ^ oo) for fixed W = {b — a). 
The barrier height for Kj < Kc is 



Gb 



27r/ioA 



In 



(I)- 



(19) 



B. Ginzburg-Landau calculation for ^ ^ W 

The calculation of the critical-current reduction factor 
R of Eq. (18) assumed that ^ <C a, but the amount of 
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the reduction is not as great when ^ « a or when ^ ^ a. 
For the latter case it is straightforward to calculate the 
critical current as an extension of the Ginzburg-Landau 
(GL) calculation^^"^° that yields the GL critical current 
density of a long, straight strip, Jgl = '/>o/3'^/^7r/xo^A^, 
critical sheet-current density, Kql = 2(/)o/3^/^7r/ioCA, or 
critical current, Iql = 2^oW^/3'^^^7r/xo^A, which assumes 
that that both d and W are much smaller than both A and 
^. For the geometry shown in Fig. 1, the GL sheet-current 
density isi^-^o ^ -(2/V/zoA)[A-h (</)o/27r)V7], where 
/ is the magnitude of the normalized order parameter. 
Since we can neglect the self- field, we may choose a gauge 
such that the current around the arc is determined by 
the gradient of the phase of the order parameter, V7 = 
—(tik/p, where fc is a dimensionless constant. The sheet- 
current density becomes K — (jyK^, where 



(20) 



Vs = po/p is the magnitude of the superfluid velocity in 

GL dimensionless imits,^^"^" and po = fc^ is a measure of 
the total current aroimd the arc. When ^ W = (b—a), 
the first Ginzburg-Landau equation yields 1 — v^, 

so that the radial dependence of K^{p) is given by 



P 



(21) 



The integral of over the width W = b — a yields 
the total current / and the average sheet-current density 
Kj = I/W, 



Ki 



( \ r Po ln(6/a) _ pl{b + a) 

Kj is maximized when po = pomax, where 

21n(6/a) 1 1/2 



POn 



= ab 



3(62 -a2) 



(22) 



(23) 



The corresponding maximum vahic of Kj is the arc's 
Ginzburg-Landau critical sheet-current density. 



KcGL = KglRgl, 
where Kql = 20o/3^/^7r/io€A, and 



Rgl = 



ab 



v/(« + &)/2 



ln(6/a)i3/2 
b — a - 



(24) 



(25) 



is the GL reduction factor due to the inhomogeneous cur- 
rent density around the arc. For (b — a) = W = a, the 
case shown in Fig. 1, Rgl = 0.94. Expanding for small 
values of (6 — a)/ (6 + a) yields Rg l ~ 1 — [(^ ~ o^)/(^ + 
a)f/2. 



VI. RADIUS-OF-CURVATURE EFFECT 

As we have seen above, when the superconducting strip 

is curved, current crowding occurs at the point of mini- 
mum radius of curvature pc-, and this reduces the critical 



sheet-current density of the strip below the value given 
in Eq. (14) when ^ ^ pc- However, for values of pc of the 
order of ^ or smaller, the critical-current reduction due to 
current crowding is ameliorated by a radius-of-curvature 
effect. In strips fabricated with sharp corners, one might 
at first expect the critical current to vanish, because the 
sheet-current density diverges at the sharp inner cor- 
ner. However, for such cases the radius-of-curvature ef- 
fect partially compensates for the current-crowding effect 
and leads to a critical current that is reduced by a factor 
proportional to (^/a)", where a is a characteristic linear 
dimension of the strip and n is a geometry-dependent 
fractional exponent. 



A. Rounded 180-degree turnaround 

To analyze this radius-of-curvature effect near the in- 
side corner of a rounded 180-degree turnaround, we use 
the conformal mapping,^^ 



C(w) = -pcw'^/2 - ipcW, 
w{0 = i(V2C/pc + 1-1), 



(26) 
(27) 



which maps points in the upper half w plane into points 

in the C plane (C, = x + iy) to the right of the parabola 
X = —y^/2pc. This parabola has radius of curvature pc 
at the origin. See Fig. 2. 

The complex potential describing uniform current flow 
from right to left in the w plane is Gw = —KoPcW, and 
in the ( plane the corresponding complex potential 



ga) = -KoPcwiO 



(28) 



describes the current crowding around the point (x, y) = 
(0,0). The imaginary part yields the stream function 
S{x,y) = 9^c(C), shown in Fig. 2(a). The complex cur- 
rent density is 



-iKn 



rfC ^/K/pc+l 



--K^{x,y)-iKy{x,y). (29) 



The constant prefactors in Eq. (28) arc chosen such that 
the current density at the origin is K{0, 0) = yKo- 

The critical current is reached when a vortex can be nu- 
cleated from the parabolic tip at {x, y) = (0, 0). To calcu- 
late it, we first need the work term Wi{ry) ~ (f)oAI{ry) = 
<j)o[S{0,0) — S{Xy,yy)], which can be obtained from Eq. 
(28). For ry = {S,0) or C,, = S, 



AI{5) = K^pcW25/p, + 1-1). 



(30) 



The self-energy Eseifivy) = (j)oIcirc{rv)/2 can be ob- 
tained as follows. The complex potential describing the 
sheet-current flowing around a vortex in the C, plane of 
Fig. 2(b) is given by Eq. (9) but with w{Q given by Eq. 
(27). We are interested only in the case when the vortex 
is on the x axis at = {5, 0) or C^^ = 5. The imaginary 
part is the stream function Sy{xy, yv',x, y) = ^QvciCv] C); 
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FIG. 3: Numerical results for (a) 5c/£, and (b) Koc (normal- 
ized to (^o/evr^oC^) as functions of the ratio of the radius of 
curvature pc (see Fig. 2) to the coherence length ^. Expan- 
sions in powers of Pc/C are shown as dashed curves for 5c [Eq. 
(34)] and Kqc [Eq. (35)]. Expansions in powers of ^/pc arc 
shown as dot-dashed curves for 5c [Eq. (36)] and Koc [Eq. 
(37)]. 



FIG. 2: (a) Current crowding around a parabolic bond of ra- 
dius of curvature pc, shown by the contour plot of the stream 
function S(x,y) = '^Qi^{x + iy) [Eq. (28)], which has the 
value S = along the parabolic boundary, x = —y^/2pc. 
The contours correspond to streamlines of the sheet-current 
density K, and the arrow shows the current direction, (b) 
Current flow generated by a vortex interacting with the 
parabolic boundary, shown by the contour plot of the stream 
function Sv(xv,y^;x,y), which has the value Sv = for 
{x, y) along the boundary. The contours, shown here for 
{xv,yv) = (Pc,0), correspond to streamlines of the vortex- 
generated sheet-current density , and the arrow shows the 
direction of the current. 



shown in Fig. 2(b). When ^ ^ a, Idrcifv) = Sv{d, 0]6 + 
^, 0) can be evaluated as described in Sec. IV. The result 
is 



-In 



2pc /25 _ / I2S 



TT/ZqA 

The Gibbs free energy is 

-2pc [26 




G 



27r/xoA 



In 




+ 1-1 



+ 1-1) 



)]• (31) 




+ 1-1 



)■ 



(32) 



The position of the barrier dj, is the value of 6 at the 



force-balance condition, dG/dS = 0. The critical sheet 
current density at {x,y) = (0,0) is reached {Kq = Kqc) 
when the barrier height G is reduced to zero at Sb = 5c, 
and 



Koc 



(33) 



where numerical results for dc/£, and fcoc as functions of 
Pc/C are shown in Fig. 3. Analytic expansions (including 
only the first few terms in the series) for Pc/C ^ 1 are 



koc 



2 



2MIJ' 



2^ J 3e^\2^r 



and corresponding expansions for Pc/^^ 1 are 



e e 

2 ^TV2p, 



\2oJ 



(2I) ' 



1 + e 



\2prJ 2 \2o.J 



e 

4 \2p. 

2 , C s2 



.2pe 



(35) 

(36) 
(37) 



Regardless of the size of pc relative to ^, vortex nucle- 
ation occurs when the barrier height is reduced to zero at 
a distance 5c of the order of ^ from the point of minimum 
radius of curvature [see Fig. 3(a)]. When ^ <C Pc, vortex 
nucleation occurs when the sheet-current density Kq at 
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FIG. 4: (a) Current crowding around a generalized hyperbolic 
bend of radius of curvature pc, shown by the contour plot of 
the stream function S{x,y) = QQ(;{x + iy) [Eq. (40)], which 
has the value 5 = along the boundary. The contours corre- 
spond to streamlines of the sheet-current density K, and the 
arrow shows the current direction, (b) Current flow gener- 
ated by a vortex interacting with the boundary, shown by the 
contour plot of the stream function 5„(a;„, y^; x, y), which has 
the value Sv = for (x, y) along the boundary. The contours, 
shown here for {xv,yv) = (Pc, 0), correspond to streamlines of 
the vortex-generated sheet-current density Ky , and the arrow 
shows the direction of the current. 



this point exceeds (po/enno^A. However, in the opposite 
hmit (pc ^ C)i vortex nucleation does not occur until Kq 
reaches much larger values [see Fig. 3(b)]. 



B. Rounded 90-degree turn 

To analyze the radius-of-curvature effect near the in- 
side corner of a rounded 90-degree turn, we use a different 
conformal mapping,^^ 

C{w) = ipMl-iwf/'-l], (38) 
w{0 = m/pc + lf/'-l], (39) 

which maps points in the upper half w plane into points 

in the C plane (C = x + iy) to the right of a generalized 
hyperbola whose radius of curvature is pc at the origin. 
See Fig. 4. 

The complex potential describing uniform current flow 
from right to left in the w plane is Gy, = —{KoPc/2)w. 



In the C plane the corresponding complex potential 

GdO = -{Kopc/2)w{0 (40) 

describes the current crowding around the point (x, y) = 
(0,0). The imaginary part yields the stream function 
S{x,y) = ?sQt^{(), shown in Fig. 4(a). The complex cur- 
rent density is 



^(0 = 



dC (3C/pc + l)i/3 



--K^{x,y)-iKy{x,y). 

(41) 

The constant prefactors in Eq. (40) are chosen such that 
the current density at the origin is If (0,0) = yKo. The 
critical current is reached when a vortex can be nucleated 
from the tip at {x,y) = (0,0). To calculate it, we first 
need the work term Wi{ry) = (/)oA/(r„) = ^o['S'(0,0) — 
S{xy,yy)], which can be obtained from Eq. (40). For 
Ty = {d, 0) or Cv = S, 



AI{S) = iKoPc/2)[{3S/p, + 1)2/3 - 1). 



(42) 



The self-energy E^eifiry) = 4>oIcirc{ry)/2 can be ob- 
tained as follows. The complex potential describing the 
sheet-current flowing around a vortex in the ( plane 
of Fig. 4(b) is given by Eq. (9), but where w{() is 
given by Eq. (39) and we are interested only in the 
case when the vortex is on the x axis at ry = (S, 0) 
or C„ = S. The imaginary part is the stream function 
Sy{xy,yy-,x, y) = '^Gy^iCvlO, shown in Fig. 4(b). When 
^ <C a, Icirc{rv) = Sy{5,0;6 + ^,0) can be evaluated as 
described in Sec. IV. The result is 



TT/ioA""!. ^ \p, 

The Gibbs free energy is 

'3S \2/3 



{f(^')"l(^')"-)]}■ 

(43) 



G = 



2npoA 

<I>oKqPc 



1/3 



(44) 



The position of the barrier 6b is the value of 5 at the 
force-balance condition, dG/dS = 0. The critical sheet 
current density at {x, y) = (0, 0) is reached {Kq = Kq^) 
when the barrier height G is reduced to zero at Si, = 5c, 
and 



(45) 



where numerical results for 6c/^ and fcoc as functions of 
Pc/C ^re shown in Fig. 5. Analytic expansions (including 
only the first few terms in the series) for pc/S, 1 are 



ei/2 /Pe\2/3 
3 6 \^) 

/Pe\-4/3 

72 \j) 



_ 

3U 
5e-3/2 

144 



(46) 
(47) 
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FIG. 5: Numerical results for (a) Sc/£, and (b) Koc (normal- 
ized to oio/eTT/io^A) as functions of the ratio of the radius of 
curvature Pc (see Fig. 4) to the coherence length ^. Expan- 
sions in powers of Pc/C are shown as dashed curves for Sc [Eq. 
(46)] and Kqc [Eq. (47)]. Expansions in powers of (,/pc are 
shown as dot-dashed curves for 5c [Eq. (48)] and Koc [Eq. 
(49)]. 



and corresponding expansions for pc/^ ^ 1 are 



Koc 



\2pJ 



Ml) 



2pc 



\2oJ 



2 \2p. 



(48) 
(49) 



Regardless of the size of pc relative to ^ , vortex miele- 
ation occurs when the barrier height is reduced to zero at 
a distance 5c of the order of ^ from the point of minimum 
radius of curvature [see Fig. 5(a)]. When ^ ^ pc, vortex 
nucleation occurs when the sheet-current density Kq at 
this point exceeds (po /eTr/^o^A. However, in the opposite 
limit (pc ^ ^), vortex nucleation docs not occur until Kq 
reaches much larger values [see Fig. 5(b)]. 



VII. 180-DEGREE TURNAROUNDS 

In this section we address the extent to which the 
critical current is reduced by a 180-degree turnaround. 
We consider the following specific examples: a sharp 
180-degree turnaround, an optimally rounded 180-degree 
turnaround, a 180-degree turnaround intermediate be- 
tween sharp and optimally rounded, a rounded 180- 
degree turnaround at the end of straight strips, and a 
sharp rectangular 180-degree turnaround. The latter ge- 
ometry was used in the experiments of Yang et al,-'^'* 
which will be discussed in more detail in Sec. XIII. 



A. Sharp 180-degree turnaround 

Consider the current flow in a strip of width a with a 
sharp 180-degree turnaround as shown in Fig. 6. Alter- 
natively, we can think of this as a wider strip of width 2a, 
cut along the x axis for x <0. The conformal mapping^^ 

d({w) 2a w 



aw TT — 1 

Ciw) = (2a/7r)(ln \/w-l + \n Vw~- 



1 



(50) 
i7r/2),(51) 

maps points in the upper half w-plane {w = u + iv) into 
the strip —a<y<am the ^-plane (C = x+iy) as shown 
in Fig. 6. The inverse mapping is glA-en by 



w{C) = ±vT^^exp(7rC7a), 



(52) 



where the upper (lower) sign holds when SjC ^ (^C > 
0). 

The complex potential is given by Eq. (7) but with 

w{() from Eq. (52) and / = Kja, and the corresponding 



complex sheet current /C^(C) 
iKy{x,y) is 



Ki 



(53) 

-00, and for < 
The streamlines 

are obtained as 



For — a < y < 0, Kx{x, y) ^ Kj as x ^ 
y < a, Kx{x,y) —Kj as x ^ — oo. 
of the sheet current K = + yKy 
contours of the stream function S{x, y) 
the imaginary part of G(;{C), shown in Fig. 6(a). 

The critical current of the 180° turnaround is reached 
when a vortex can be nucleated at the sharp point in 
the center of Fig. 6(a) or (b) at {x,y) = (0,0), which 
corresponds to C = and w = 0. To calculate the critical 
current, we examine the behavior when a vortex is at 
C„ = S, where (5 <C a. Expanding Eq. (50), we obtain 
^^(Cw) = «(7r5/a)^/^ to lowest order, and from Eqs. (7), 
(10), and (11) we obtain the Gibbs free energy, 

-)-2^oi^,(-) . 



G = 



2irfj,oA 



In 



(54) 



Following the steps that led to Eq. (14), we obtain with 
Ki = Kc, 



^1/2 ^ 90 

^ 2p,oAKi,/iw, 



, Sc = e^il^ = 1.85^, (55) 

7re\i/2 



-i?, where R 



(?) 



(56) 



eTT/io^A 

is the reduction factor due to current crowding near 
{x,y) = (0,0). For example, R = 0.56 when ^ = 0.1a 
and R = 0.18 when ^ = O.Ola. The barrier height for 
Kj < Kc is 



Gb = 



ttpqA 



(57) 



Note that the prefactor is larger than that in Eq. (15) by 
a factor of 2, which arises from the term proportional to 
(5^/2 in Eq. (54). 
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FIG. 6; (a) Current flow in a strip carrying current Kia 
around a 180-degree turn, shown by the contour plot of the 
stream function S(x,y) = ^Qq{x + iy), which has the val- 
ues 5 = along the lines y = ±a at the outer boundaries 
and S = Kia on either side of the narrow gap along the line 
y = Q for a; < 0. The contours correspond to streamlines 
of the sheet-current density K, and the arrows show the di- 
rection of the current. The dashed curve, which corresponds 
to S = Kia/2, separates the current-crowding region close 
to {x, y) = (0, 0) from the current-expanding region outside. 
The magnitude of K is constant [K — Ki) along the dashed 
curve, (b) Vortex-generated current flow, shown by the con- 
tour plot of the stream function Svix^, yv',x, y), which has the 
values Sv = for (x, y) along the boundaries. The contours, 
shown here for {x^,y„) = (0.1a, 0), correspond to streamlines 
of the vortex-generated sheet-current density , and the ar- 
row shows the direction of the current. 



B. Optimally rounded 180-degree turnaround 

For the case of the 180° turn shown in Fig. 6, the 
complex sheet-current density is given by Eq. (53). Ex- 
amination of Fig. 6(a) reveals that current crowding oc- 
curs, i.e., K — \K\ increases along streamlines inside 
the dashed curve as the current turns around the end 
of the gap, reaching a maximum at {x,y) = (0,0). On 
the other hand, K decreases along streamlines outside 
the dashed curve in the current-turnaround region. The 
dashed curve, which corresponds to the contour for which 



S = Kja/2, but which also can be obtained by setting 
K = Kj, is given by either of the following equations: 

Xo{y) = (a/7r)ln[2cos(7r?//a)], (58) 
yo{x) — ±(a/7r) cos^"^[exp(7rx/a)/2]. (59) 

In the latter equation, x < (a/7r)ln2 — 0.221, and the 
upper (lower) sign holds for positive (negative) values of 
y. Note that y in Eq. (59) rapidly approaches ±a/2 for 
X < —a as exp(7ra;/a) — > 0. 

The above results can tell us the optimal film design 
that will prevent any significant reduction of the criti- 
cal current due to current crowding at a 180-degree turn. 
Consider a long superconducting strip of width W and 
critical sheet current given by Eq. (14). If we wish the 
direction of the current to change by 180° at the corner 
of a strip for which the outer boundaries are straight, as 
shown in Fig. 6(a), the inner boundary of the strip should 
be chosen to be the smooth curve given by Eq. (58) or 
(59) but with a/2 = W. The minimum radius of curva- 
ture of this curve is pc = 2W/tt = 0.637VF, which occurs 
at {x,y) = (a;o(0),0) = (2VK In 2/7r, 0) (0.44W,0). 
As discussed in Sees. V and VI, so long as ^ ^ pc, the 
self-energy Egei / of a nucleating vortex for small 6 is the 
same as in Eq. (12) to excellent approximation. More- 
over, since the sheet-current density along the entire in- 
ner boundary is constant with the value Kj, the work 
term Wi and hence the entire Gibbs free energy for small 
d are very nearly the same as in Eq. (12). Therefore, the 
critical current for a strip with a 180-degree turnaround 
of the above-described design, i.e., the area between the 
dashed curve and the outer boundary in Fig. 6, should be 
the same as that of a very long strip of constant width W 
[see Eq. (14)], so long a.s £, <^W. See also Sec. VII D 3, 
where the optimally rounded 180-degree turnaround is 
examined from a different perspective. 

C. 180-degree turnaround intermediate between 
sharp and optimally rounded 

Consider a long strip of width W whose critical sheet 
current is given by Eq. (14) when ^ <^ W. In Sec. VII A, 
we showed that current crowding reduces this critical 
current by the factor R = (tt^/W^)-^/^ at the sharp 180- 
degree turn as shown in Fig. 6 when a = W. In Sec. 
VII B, we described an optimal 180-degree turnaround 
geometry that avoids current crowding, such that there 
is no reduction of the critical sheet current given by Eq. 
(14) so long as ^ ^ W. However, the optimal 180-degree 
turnaround requires that the gap between strips be dou- 
ble the strip width W, so that in meander arrays the 
filling factor (fraction of surface area covered by the su- 
perconducting film) is only 1/3. Device designers may 
wish to increase the filling factor at the price of reducing 
the critical current. In this section we therefore show how 
to estimate the critical current in 180-degree turnarounds 
that are intermediate between the sharp and optimally 
rounded cases discussed in Sees. VII A and VII B. 
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FIG. 7: Critical-current reduction factor R [Eq. (63)] vs 
g/2a = (a — W)/a for several values of f/o. 



Consider a 180-degree turnaround consisting of a su- 
perconducting strip whose shape is chosen such that 
the outer boundaries are the straight lines at y — ±a 
as shown in Fig. 6(a) and the inner boundary is the 
curve defined by the stream-function contour S{x,y) = 
^Gdx + iy) = KjW. Far to the left of the turnaround, 
the film has a nearly constant width W extending from 
y = —a to y = —{a — W) for y < 0, where it car- 
ries a sheet-current density K = xKj, and it has a 
nearly constant width W extending from y = +a to 
y = +{a — W) for y < 0, where it carries a sheet- 
current density K = —xKj. Analysis of the contour 
S{x,y) = KjW reveals that it intersects the x axis at 
{x,y) = {xw,Q), where 



xw = (2a/7r)ln[l/sin(7rW/2a)], 



(60) 



the radius of curvature is 



p^{xw) = — (1 - e-""'^/") = — cos^{nW/2a), (61) 

TT TT 

and the sheet-current density is, from Eq. (53), 

0) = ^ e^W«-l = tan(7rW^/2a). (62) 

We now make use of the results of Sec. VI, in which 
the critical sheet-current density Kqc at the point of 
minimum radius of curvature pc is given by Eq. (33), 
where /cqc is the function of Pc/^ plotted in Fig. 3(b). 
Since Kj = at the critical current, we have Kqc = 
ii'ctan(7rW/2a) from Eq. (62), such that 



K, 



R, where R ■ 



(63) 



and the argument of fcoc [see Fig. 3(b)] is pc/S, = 
(2a/7r^) cos^(7rH^/2a), obtained from Eq. (61). Figure 
7 shows calculated values of i? as a function of g/2a, the 
ratio of the width of the gap g = 2a — 2W to the distance 
2a between the two outer boundaries in Fig. 6. R de- 
creases monotonically as g decreases and reduces to Eq. 
(56) in the hmit as 5 — )■ and W ^ a. 



As expected, for very small values of ^/a (see plot for 
$/a = 0.001), R = 1 when g/2a 1/2 or = a/2. For 
larger values of ^/a, on the other hand, the behavior of 
Rvs g/2a becomes more strongly affected by the radius- 
of-curvature effect (see Sec. VI), which counteracts the 
current-crowding reduction of Kc- The portions of the 
curves for which R > 1 tell us the values of ^/a, g, and 
W for which the critical current is limited by vortex nu- 
cleation not at the point of minimum radius of curvature, 
{x, y) = {xw, 0), but rather somewhere along the straight 
portion of the strip, far from the bend. 



D. Rounded 180-degree turnaround at the end of 
straight strips 

We now use a method inspired by that of Cockroft^ to 
calculate the current-crowding critical-current reduction 
factor R for the case that straight strips with constant 
width W and gap 2r' between them arc connected at their 
ends with a rounded corner, as shown in Fig. 8. The re- 
sults depend in detail upon the geometry chosen. In Sec. 
VII D 1 we discuss the behavior for which the dimensions 
shown in Fig. 8(a) obey W = W, in Sec. VII D 2 we 
examine the behavior for W 00 and r' /W < 1, and 
in Sec. VII D 3 we treat the limiting case for W' 00 
and r' /W — 1, which produces a curving inner bound- 
ary of the turnaround corresponding to optimal rounding 
discussed in Sec. VII B. 



1. General case 
The conformal mapping^ 

dCiw) 



aw) 



dw 

-A 



A——^=—=-, (64) 
— WW — ay/w — 



2i 



--[{a - l)n((/)i,ni,fci) 



V(a+1)(6-1) 

-{a + l)F{(bi,ki)] 
\/a — w + \fb 



j b+l f y'a-w + ^/b-w \\ 



W/g{b), 




a + V 



a + V 



h = 



' 2(6 - a) 
{a+l){b-iy 



(65) 
(66) 
(67) 

(68) 

(69) 

(70) 



maps points in the upper half w-planc {w = u + iv) into 
the area ABCDEF in the upper half ^-plane (^ = a; -|- iy) 
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shown in Fig. 8. The inverse mapping can be ob- 
tained numericaUy. Here, F((j),n,k) and Il{(f>,n,k) are 
elliptic integrals of the first and third kind with param- 
eter n and modulus fc,^^ and K{k) = F{'K/2,k) and 
n(n, k) = n(7r/2, n, k) are the corresponding complete el- 
liptic integrals. Special points in the ^ and w planes (see 
Fig. 8) are related as follows: A, w = — oo; B, w = —1; 
C, w = +1; D, w = a; E, w = b; and F, w = +oo, where 
1 < a < 6. The lengths CD = W' , r, and r' shown in Fig. 
8(a) are determined as functions of a and b as follows: 



W 



= ^{ J.,,, \{b+l)K(k,) 
-(6-a)n(n3,fc2)] 



+2 



b-l 



In I 



/IT 



)}• 



b-l 



(a+l)(6-l)' 



= 2A 



b-l 
2A 



m 



V(a + 1)(6-1) 



1 V6 - 1 

[{l + a)K{ki) 



(71) 
(72) 

(73) 



+(l-a)n(ni,A;i)]. (74) 

Note that DE = VF + r'. 

The complex potential describing the current flow 
within the area ABCDEF shown in Fig. 8 is ^^(C) = 
Qw{w), where 




x/w 



g^{w) = - ln{2w - 1 - a + 2y/w - ly/w - a) (75) 

TT 

and w = w{(). The imaginary part of this is the 

stream function S, whose contours, shown in Fig. 8, 
are the streamlines of the sheet-current density K = 
xKx + yKy. The corresponding complex sheet current, 
/C^(C) = dg^{Q/dC, = K^{x,y) - iKy{x,y), is 



^c(C) 



lA 



w ■ 



(76) 



For r' < y < W + r' and x ^ —W, which corresponds to 
w oo, we have to good approximation = —Kj = 
—lA/gib) = —I/W. However, around the arc BC, which 
corresponds to w = u, where — 1 < u < 1, we find that 
-K'bc = |/C| is given by 



Kbc ^ ^/b+l + y/b^ 
Ki ~ ^/2 



(77) 



Since b is required to obey 6 > 1 except in limiting cases, 
this equation tells us that current-crowding (i^ec > Kj) 
almost always occurs. As discussed in Sees. V and VI, so 



FIG. 8: (a) Current flow in a film carrying current KjW 
from a straight strip of width W around a rounded 180-degroo 
turn of length CD = W' = W into another straight strip of 
width W, shown by the contour plot of the stream function 
S{x, y) = QG({x + iy), which has the values 5 = along the 
outer boundaries of the film and S = KjW along the inner 
boundary. Here r' — 0.50014^ and r = 0.402VF. The contours 
correspond to streamlines of the sheet-current density K, and 
the arrows show the direction of the current, (b) Same as (a), 
except that the connection is infinite in length (CD = oo), 
r' = O.SOOW^, and r = 0.525^^. 



long as ^ ^ r', this current-crowding therefore leads to 
a critical sheet-current density of the form 



-R, where R ■ 



V2 



Vb+1 + Vb^ 



(78) 



is the reduction factor due to current crowding at the 
curving inner boundary BC, relative to the K^. for a long 
straight strip [Eq. (14)]. However, to treat the case for 
which r' < an approach similar to that in Sec. VI 
would be required. 

Using Eqs. (66), (71), and (74), one can obtain two 
equations from which the values of b and a can be ob- 
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tained for desired values of the ratios W /W and r' /W. 
Figure 8(a) was calculated for the following parameters: 
W = W = 1, r' = 0.500, r = 0.402, A = 0.102, 
a = 1.216, and b = 1.572. 

Combining displaced mirror images of turnarounds like 
these in a two-dimensional layout (a "boustrophedonic" 
pattern) results in a fill factor / — W/p, where p = 
W + 2r' is the pattern period (pitch) and W and 2r' 
are the strip and gap widths away from the turnarounds. 
For example. Fig. 8 corresponds to / = 1/2. The solid 
curves of Fig. 9 show plots of R, r' /W = (1 — f)/'2,f, 
r'/r, and 6-1 vs / = l/(l + 2r7VF) for W /W = 1. For 
/ <C 1, the following expansions have been used to plot 
the functions for / < 0.12: 



2.0 



5b 
R 



b-l = 



128 , 7r(l-/), 
— exp ^ ^1 



5bl2 



2/ 

,/5bj2), 



1-/ 



1 - /[I + (4/7r) ln(4/7r)] ' 



(79) 
(80) 

(81) 



While Fig. 9 shows R approaching zero for / ^ 1 and 
r'/W 0, bear in mind that Eq. (78) is valid only for 



0.5 



no 



r'/W \ 

r'/r 


i/=l/3 


/=i/2 / 

b-l jb-l 
/ ' 









0.2 



0.4 



0.6 



0.8 



1.0 



/ 



FIG. 9: Plots of the current-crowding critical-current reduc- 
tion factor R, 6b = b—1, and the ratios r'/W and r'/r vs the 
fill factor / = !/(! + 2r' /W) for W = W (solid curves) and 
W' = oo (dashed curves). The vertical dotted line marks the 
parameters for Fig. 8(a) (/ = 0.5, r'/W = 0.5, R = 0.599, 
6 - 1 = 0.572, and r'/r =1.243) and Fig. 8(b) (/ = 0.5, 
r'/W = 0.5, R = 0.707, 6 - 1 = 0.250, and r'/r =0.953). The 
vertical dot-dashed line marks the smallest / for which solu- 
tions can be found for W' = oo and the value of / = 1/3 for 
which r'/W = 1, r/W — > oo, and the inner boundary takes 
the optimally rounded shape for a 180-degree turnaround so 
that R=l. 



2. W' ^oo and r'/W < 1 

In the limit CD = W' ^ oo,b ^ a, and Eqs. (64)-(74) 
simplify to 

C_'(w) = = -A , Y; ^ ,(82) 



dw 



\/w — l{w — a) 

1 - ?i; la + 1 



A = W/g{a), 

a + 1~ 

a — 1 - 

^Inf^), 

o-i y^^^J' 

a+l 



g{a) = 7r(l 
r = 2A 



(83) 
(84) 
(85) 

(86) 

(87) 



For W — >■ 00, Eqs. (84) and (87) can be used to obtain 
an equation from which b = a can be determined as a 
function of r'/W, but only for r'/W < 1. Equations 
(75)- (78) still apply, and Fig. 8 (b) was calculated for the 
following parameters: W = 1, W' — >■ oo, r' = 0.500, r = 
0.525, A = l/47r, and b = a = 1.25. The dashed curves 
of Fig. 9 show plots of R, r'/W = (1 - f)/2f, r'/r, and 
6-l = a-lvs/ = l/{l + 2r'/W) for W' /W = 1. While 



the dashed curve for R in Fig. 9 shows R approaching zero 
for / -J> 1 and r'/W -)> 0, recaU that Eq. (78) is vahd 
only for r' ^ ^. 



3. Optimal rounding when W' oo and r/W oo 

Note from the dashed curves in Fig. 9 that in the limit 
W' -> oo, we find R = 1, r'/W = 1, r/W -> oo, and 
/ = 1/3 in the limit a = 6 — >■ 1. The dependence of 
({w) then becomes dominated by values of w very close 
to 1 and a, and it is appropriate to introduce the variable 
to = {w-l)/{a-l) in Eqs. (82)-(87) and to take the limit 
as a ^ 1. This leads to the conformal mapping 



C'('^) 
C('^) 

wiC) 



dC{oj) __W{l + ^/cJ) 

doj \/oj{io — 1) ' 

[1 - exp(-7rC/2VK)]2, 



(88) 

(89) 

(90) 



which maps points in the upper half w-plane into the area 
BCDEF in the upper half C-plane {( ~ x + iy) shown in 
Fig. 10. Special points in the C and uj planes (see Fig. 
10) are related as follows: B, w = — oo; C, w = 0; D, 
w = 1 — e; E, w = 1 + e; and F, w = +oo, where e is 
a positive infinitessimal. The curves BC and B'C for 
a; < arc given by 

yopt{x) = ±(2M^/7r) cos-i[exp(7rx/2W')] (91) 
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FIG. 10: Current flow around a 180-degree turnaround with 
an optimally rounded inner boundary, shown by the contour 
plot of the stream function Six, y) — QQ({x + iy), which has 
the values S = G along the outer boundaries of the film and 
S = KiW along the inner boundary. The contours corre- 
spond to streamlines of the sheet-current density K, and the 
arrows show the direction of the current. Since there is no cur- 
rent crowding along the curve BCC'B', the critical current is 
predicted to be the same as for a long, straight strip. 



with the upper (lower) sign holding for BC (B'C). These 
curves correspond to the optimally rounded inner bound- 
ary of a 180-degree turnaround discussed in Sec. VII B. 
Note that yopt{x) in Eq. (91) is the same as yo{x) in Eq. 
(59) with a = 2W and the origin shifted along the x axis 
by Ax = (a/7r)ln2. 

The complex potential describing the current flow 
within the area BCDEF shown in Fig. 10 is 

g^c) = - ln[^^;(c) + ^^(c^^)]. (92) 

TT 

The imaginary part of this is the stream function S, 
whose contours, shown in Fig. 10, are the streamlines 

of the sheet-current density K = xK^ + yKy. The cor- 
responding complex sheet current, /C^(C) = dGdO/'^C = 
Kx{x,y) - iKy{x,y), is 



W^l + V^' 



(93) 



For yopt < y < and x <C —W, which corresponds 
to oj{C) oo, we obtain = —Kj, where Kj = I/W. 
Along the arc BC, which corresponds to w = u, where 
—00 < u < 0, we find that Kbg = |^| = Kj. In other 
words, there is no current crowding, and the critic;al cur- 
rent is predicted to be exactly the same as for a long 
straight strip [Eq. (14)]. As discussed above, however, 
this conclusion is based on the assumption that the co- 
herence length ^ is much smaller than the radius of cur- 
vature at the origin pc = 2W/'k. To treat the case for 
which ^ > an approach similar to that in Sec. VI 
would be required. 
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FIG. 11: (a) Current flow in a film carrying current KiW 
around a sharp rectangular 180-degree turn, shown by the 
contour plot of the stream function S{x,y) = QGci^ + w)i 
which has the values S = along the lines y = ±a at the 
outer boundary of the film and S = KjW along the inner 
boundary [\y\ = ft = a — for a; < or [yl < /i for x = 0]. 
(In this figure, h ^ a/2 = W, and a ^ V3/2 = 0.866.) 
The contours correspond to streamlines of the sheet-current 
density K, and the arrows show the direction of the current, 
(b) Vortex-generated current flow, shown by the contour plot 
of the stream function Sv{xv,yv; x,y), which has the values 
3.^=0 for (x, y) along the boundaries. The contours, shown 
here for {xv,yv) = (0.05a, —0.55a), correspond to streamlines 
of the vortex-generated sheet-current density Ky, and the 
arrow shows the direction of the current. 



E. Sharp rectangular 180-degree turnaround 

We next examine the current flow around a sharp rect- 
angular turnaround shown in Fig. 11. The fllm (width 
2a) is centered on the x axis, but a slot (width 2h) with 
sharp 90-degree corners has been cut out of the center for 
a; < 0. The lower arm of width W = a — h carries a uni- 
form sheet-current density K = xKj until it reaches the 
turnaround. The current then turns around and finally 
fiows in the upper arm with current density K = —xKj. 
We wish to calculate the critical current at which the first 
vortex is nucleated at one of the sharp inner corners of 
the turnaround. 
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The conformal niapping^° 

d({w) 2a s/w^ — o? 



C{w) 



dw 
'2a\ 



CM = (^)[cosh-(^)-^./2 



. _i /w^(2 - a^) - 



cosh" 



-1 (un^ 



- 1) 



(94) 



,(95) 



where a = yj2{h/a) — (h/a)"^, maps points in the upper 
half w-plane {w = u + iv) into the area —a<y<a for 
X > OT h < \y\ < a for a; < in the C-plane = x + iy) 
as shown in Fig. 11. The inverse mapping w(C) can be 
obtained numericahy. 

The complex potential GciC) is given by Eq. (7) using 
the inverse mapping w{Q with I = KjW, and the cor- 
responding complex sheet current /C^(C) = dG,^{()/d( = 
Kx{x,y) - iKy{x,y) is 



^c(C) 



(96) 



For —a <y < —h, Kx{x,y) Kj as x ^ — oo, and for 
h < y < a, Kx{x,y) — > —Kj as x — > — oo. The stream- 
lines of the sheet current K = + yKy are obtained as 
contours of the stream function S{x,y) = 'QQ^ix + iy), 
the imaginary part of GciC): shown in Fig. 11(a). 

The critical current of the rectangular 180° turnaround 
is reached when a vortex can be nucleated at one of the 
sharp inner corners shown in Fig. 11(a) or (b). The lower 
right corner at C = —ih corresponds to w = a. Ex- 
panding Eq. (50) for = {S/V2, h - 6/^2) or = 
—ih + S/y/2 — iS/y/2 along the diagonal extending away 
from the lower right inner corner of Fig. 11(a), we obtain 



wiCv) 



/W'xVa /37r(5\2/3 



(2a)V3 \ a J V 4o 7 



(97) 



to lowest order, where 5 <C a. From Eqs. (7), (10), and 
(11) we obtain the Gibbs free energy, 



\iTaJ 



(98) 



Following the steps that led to Eq. (14), we obtain with 



~47r2/3/XoAii:/^V3 



, 5,= — A = 1.49^, (99) 



K.= -^R, where R = -f ^)'^' (100) 

is the reduction factor due to current crowding at one 

of the sharp inner corners, and (because h = a — W) 
a = -^/l — {W/a)^. The barrier height for Kj < Kc is 



G.= ^ln^^' 



47r/xoA 



(101) 



Note that the prcfactor is larger than that in Eq. (15) by 
a factor of 3/2, which arises from the term proportional 
to ^2/3 in Eq. (98). 

In the limit when ^ a, a — >■ 1, in which case Eq. 
(100) yields R « {3/2){n^/2Wy/^, a result nearly the 
same as that in Eq. (108) with a = W but larger by a 
factor of 2^/^ = 1.26 because of the geometry differences. 

The steps leading to Eq. (100) should be valid so long 
as h = {a — W) <^ ^, but the above approximations fail 
for very small gap widths q = 2h ^ ^. In this case, 
a ~ i^/Wy/'^, and Eq. (100) yields R ~ {^/Wy/^. In 
the limit as /i ^ 0, the value of R is given by Eq. (56). 



VIII. 90-DEGREE TURNS 

In this section we first calculate how much the criti- 
cal current is reduced by a sharp 90-degree turnaround, 
and we then describe the shape of an optimally rounded 
90-degree turn that should exhibit no critical-current re- 
duction. 



A. Sharp 90° turn 

Consider the current flow in a strip of width a with 
a right-angle turn as shown in Fig. 12. The complex 
potential that describes the flow of current is 



^;c(C) = -— inHO-i], 

TT 



(102) 



where the conformal mapping from the w plane {w 
u + iv) to the C, plane [Q = x + iy) is^^ 



Q'{w) = 



dC{w) aVT 



w 



(103) 



dw 7r(l — w)^/w^ 
C(w) = -[cos-i(l-M;)-Fcosh-i f^— )]/104) 

IT I \ 1 — J 

and w(C) is the inverse of C{w), which is readily obtained 
numerically. The imaginary part of G is the stream func- 
tion, S{x, y) = "^Gix + iy). Shown in Fig. 12(a) is a con- 
tour plot of S{x, y) whose contours correspond to stream- 
lines of the sheet-current density K. 

To calculate the critical current at which a vortex is 
nucleated from the inner corner, note that the corner at 

= a + ia corresponds to w = 2. Expanding Eq. (103) 
about this point with (^y = a + ia + 6e~^'"^^ yields 



w{Cy) = 2- i(37r(5/aV2)2/^ 



(105) 



We can use this expansion following the procedure of Sec. 
IV to obtain the Gibbs free energy. 



G = 



27r/xoA 



In 



(f)-^oi^.(|^)"^(3^)^/^. (106) 
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Ki < Kc is 
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FIG. 12: (a) Current flow in a strip carrying current Kia 
around a 90° turn, shown by the contour plot of the stream 
function ^(a;, y) — QQ{x+iy) [Eq. (102)], which has the values 
S — along the inner boundary {y = a for x > a and x — a 
for y > a) and S = —Kja along the outer boundary {y — 
for X > and x = for y > 0). The contours correspond 
to streamlines of the sheet-current density K, and the arrows 
show the direction of the current. The dashed curve, which 
corresponds to 5* = —Kia/2, separates the current-crowding 
region near the inner corner from the current-expanding re- 
gion near the outer corner. The magnitude of K is constant 
{K = Ki) along the dashed curve, (b) Vortex-generated cur- 
rent flow, shown by the contour plot of the stream function 
Sv{xv,yv',x,y), shown here for {xv,yv) ~ (0.9a, 0.9a), which 
has the values — for {x, y) along the inner and outer 
boundaries. The contours correspond to streamlines of the 
vortex-generated sheet-current density Kv 



Following the steps that led to Eq. (14), we obtain with 
Ki = K,, 



37r\i/3 



00 



'0 /';37r\ 

AKi\Aa) 



3 /7r$\i/3 
^ . R, where i? = - — 
eTT/io^A 2V4a/ 



1.49CC107) 



(108) 



is the reduction factor due to current crowding at the 
sharp inner corner. For example, R — 0.64 when ^ = 0.1a 
and R = 0.30 when ^ = O.Ola. The barrier height for 



Gb 



In 



(109) 



Note that the prefactor is larger than that in Eq. (15) by 
a factor of 3/2, which arises from the term proportional 
to (5^/3 in Eq. (106). 



B. Optimally rounded 90° turn 

For the case of the right-angle turn shown in Fig. 12(a), 
the complex sheet-current density in the C plane obtained 
from Eq. (102) via K^{C) = dg^{C)/d( is 



iKy = 



(110) 



where K, the magnitude of JC(^, is also the magnitude of 

K; i.e., K = |/Ccl = \K\ = ^K^+K^, and K varies 

along all of the contours of Fig. 12(a) except one. For 
contours close to the inner corner, K increases upon ap- 
proaching the corner, which we call current crowding; for 
the contour along the inner boundary, K even diverges at 
the corner. As shown in Sec. VIII A, this current crowd- 
ing leads to a significant reduction in the critical current, 
because vortices nucleate preferentially at the sharp in- 
side corner. On the other hand, for contours close to the 
outer corner, K decreases upon approaching the corner, 
which we call current expansion. 

However, as first discovered by Hagedorn and Hall,^ 
there is one special contour near the middle for which K 
remains constant with the value Kj; this contour, which 
we call the optimal contour, is shown as the dashed curve 
in Fig. 12(a). From Eq. (110) we find that in the w 
plane the optimal contour corresponds to w = 1 + iv, 
where v > 0, such that the optimal contour in the C 
plane, Co{v) = Xo{v) + iydv), is given by the parametric 
equation 



1 H — sinh 
vr 







0] 


+ i 







1 



■ sinh ^ V 



}. (Ill) 



Alternatively, v can be eliminated to obtain as a func- 
tion of a; > a/2 or Xo as a function of y > a/2: 



Voix) = "|l + -sinh ^ 



r 2 _i 

<^ 1 + - sinh 

L TT 



1 

sinh[(7r/a)(x — a/2) 
1 



sinh[(7r/a)(y — a/2) 



};i2) 
}ll3) 



The above results have important consequences, be- 
cause they can tell us the optimal film design that will 
prevent any significant reduction of the critical current 
due to current crowding at a right-angle turn. Consider 
a long superconducting strip of width W and critical 
sheet current given by Eq. (14). If we wish the direc- 
tion of the current to change by 90° at the corner of 
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a strip for which the outer boundary is a right angle, 
as shown in Fig. 12(a), the inner boundary of the strip 
should be chosen to be the smooth curve given by Eqs. 
(111)-(113) but with a/2 = W. The minimum radius of 
curvature of this curve is AW/n = 1.27W, which occurs 
at Xo = Vo^ W[l + (2/7r)sinh-^(l)] = 1.56W. As dis- 
cussed in Sees. V and VI, so long as ^ is much smaller 
than the minimum radius of curvature, the self-energy 
Eseif of a nucleating vortex for small 5 is, to excellent 
approximation, the same as in Eq. (12). Moreover, since 
the sheet-current density along the entire inner bound- 
ary is constant with the value if/, the work term Wj and 
hence the entire Gibbs free energy for small 6 arc very 
nearly the same as in Eq. (12). Therefore, the critical 
current for the strip with a corner of the above-described 
design, i.e., the area between the dashed curve and the 
outer boundary in Fig. 12(a), should be the same as that 
of a very long strip of constant width W [see Eq. (14)], 
so long as ^ <C W. 



C. Rounded 90-degree turnciround at the end of 
straight strips 

We now use Cockroft's method^ to calculate the 
current-crowding critical-current reduction factor R for 
the case that two straight strips with constant width W 
are connected by a 90-degree rounded corner, as shown 
in Fig. 13(a). The results depend upon the parameters 
chosen. In Sec. VIII C 1 we discuss the behavior for which 
the "radius" r shown in Fig. 13(a) is finite, and in Sec. 
Vine 2 we treat the limiting case for r/W ~> oo, which 
produces a curving inner boundary corresponding to op- 
timal rounding discussed in Sec. VIII B. 



1. General case 
The conformal mapping^ 

dC(w) Vw + i + y^Vw' - 1 



A 
Ciw) 



dw {w - a)\/w - b 

W/g, a={b^ + l)/2b, 



(114) 
(115) 



2A 



In 



lb+1 1 



'^/w + 1 + \/w - b 
' TT+T 

~ \b+l){w-b) 



-htan 



6-1 y(6-l)(w-hl) 
\/w — 1 + — b\ 

1 Ub-l){w-b)- 



{b+l){w-l) 



(116) 
(117) 
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FIG. 13: (a) Current flow calculated in Sec. VIII C 1 for a 
film carrying current KjW from a straight strip of width W 
around a rounded 90-dcgree turn into another straight strip 
of width W, shown by the contour plot of the stream func- 
tion S{x, y) = QQ(;{x + iy), which has the values S = along 
the outer boundaries of the film and S = KjW along the 
inner boundary. The contours correspond to streamlines of 
the sheet-current density K, and the arrows show the direc- 
tion of the current. Here r/W = 1/2, and current crowding 
along the arc BC results in a critical current reduction fac- 
tor R = 0.654. (b) Current flow calculated in Sec. Vlll C 2, 
where the connection is infinite in length and r/W — ^ oo, a 
geometry for which there is no current crowding along the arc 
ABC and no critical-current reduction. 



maps points in the upper half w-plane {w = u + iv) into 
the area ABCDEFG in the upper half C-plane [Q = x + 

iy) shown in Fig. 13(a). The inverse mapping w(C) can 
be obtained numerically. Special points in the C, and w 
planes are related as follows: A, it; = — oo; B, w = —1; C, 
w = +1; T),w = a— e; E, w = a+e; Y,w = b\ and G, w = 
-hoo; where 1 < a < 6 and e is a positive infinitessimal. 
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The "radius" r shown in Fig. 8(a) is determined from 



6+1^^^1 + ^2/(6+1) 



6-1 



VVib+1) 



) 



—2 tan 



1 



(118) 



The complex potential describing the current flow 
within the area ABCDEFG shown in Fig. 13 is 



gdO = -Hw{0-a]. 



(119) 



The imaginary part of this is the stream function S, 
whose contours, shown in Fig. 13, are the streamlines 
of the sheet-current density K = xK^ + yKy. The cor- 
responding complex sheet current, /C^(C) = dQ,^{Q/d( = 
Kx{x,y) - iKy{x,y), is 



/Cc(C) = Ki 



{Vb + T + Vb^)V^' 



^{b + i){w-i) + ^{b-i){w + iy 



(120) 



where Kj = I/W. For Q < y < W and x » W, which 
corresponds to — > oo, wc have to good approxima- 
tion Kx = Kj. However, around the arc BC, which cor- 
responds to w = u, where — 1 < u < 1, we find that 
Kbg = |/C| is given by 



Kbc Vb+1 + Vb-l 



Ki 



%/2 



(121) 



Since 6 > 1 except in the special limit 6 — >■ 1, this equa- 
tion tells us that current-crowding {K-qc > Kj) almost 
always occurs. As discussed in Sees. V and VI, so long as 
^ 4; r, this current-crowding therefore leads to a critical 
sheet-current density of the form 



K, 



R, where R 



V2 



Vb + l + Vb^ 



(122) 



is the reduction factor due to current crowding at the 
curving inner boundary BC, relative to the Kc for a long 
straight strip [Eq. (14)]. However, to treat the case for 
which r' < an approach similar to that in Sec. VI 
would be required. 

Figure 14 shows plots of R, r/W, 6—1, and a— 1 vs 6—1. 
While the plot shows R approaching zero for r/W — > 0, 
bear in mind that Eq. (122) is valid only for r ^ ^. Note 
that to achieve the goal of minimizing current crowding 
and maximizing the critical current, one should choose 
the shape of the 90-degree turn make r/W as large as 
possible. To achieve R ~ 1 requires r/W ^ oo and 
6 — 1. We discuss this limit in the following section. 



2. Optimal rounding when r/W — >■ oo and & — >■ 1 

When r/W oo and 6 — > 1, the dependence of C(w) 
becomes dominated by values of w very close to 1 and 




FIG. 14: Plots of the current-crowding critical-current reduc- 
tion factor R [Eq. (122)], r/W [Eq. (118)], 6 - 1, and a - 1 
vs 6 — 1 for the 90-degree turn discussed in Sec. VIII CI. 
The vertical dotted line marks the parameters for Fig. 13(a) 
{R = 0.654, r/W = 0.5, b = 1.383, and o = 1.053). 



6, and it is appropriate to introduce the variable w = 
{w - l)/(6- 1) in Eqs. (114)-(118) and to take the limit 
as 6 1. This leads to the conformal mapping 



C'('^) 

C(^) 



dCH ^ w{i + y/uj) 

dio oj\/oj — 1 ' 
+ cosh" 



2Wf _i 

cos 

TT V 



/UJ 



(123) 

(124) 



which maps points in the upper half w-plane into the area 
ABCDEF in the upper half ^-plane (C = a; + iy) shown 
in Fig. 13(b). Special points in the ( and w planes are 
related as follows: A, co = — oo; B, co = —1; C, oj = — e; 
D, a; = +e; E, w = 1; and F, w = +oo, where e is a 
positive infinitcssimal. The inverse mapping a;(C) can be 
obtained numerically. The inner boundary curve ABC 
for w < is given by Eqs. (111)-(113) with v replaced by 
l/-\/|ajf and a replaced by 2W . This curve corresponds 
exactly to the optimally rounded inner boundary of the 
90-degree turn discussed in Sec. VIII B. 

The complex potential describing the current flow 
within the area ABCDEF shown in Fig. 13(b) is 



^aC) = -lnu;(C)- 



(125) 



The imaginary part of this is the stream function S, 

whose contours, shown in Fig. 13(b), are the streamlines 
of the sheet-current density K = xK^ + yKy. The cor- 
responding complex sheet current, /C^(C) = dg,^{Q/d( = 
Kx{x,y) - iKy{x,y), is 



^c(C) = 



/ v/^(C) - 1 

w^i + v^' 



(126) 



For y < W and x ^ W, which corresponds to a;(C) 
±00, we obtain = Kj, where Kj = I/W. Along the 
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arc ABC, which corresponds to a; = u, where — oo < u < 
0, we find that Kbc = |^| = Kj. In other words, there is 
no current crowding, and the critical current is predicted 
to be exactly the same as for a long straight strip [Eq. 
(14)]. As discussed above, however, this conclusion is 
based on the assumption that the coherence length ^ is 
much smaller than the radius of curvature at the origin 
Pc = 2W/-K. To treat the case for which ^ > pc, an 
approach similar to that in Sec. VI would be required. 



IX. T INTERSECTION 

To measure the critical current of a narrow strip, it is 

often the case that the current is fed in and taken out 
using wide contact pads at the ends of the strip, and the 
voltage is measured using sidebar contacts that intersect 
the strip in T intersections. We defer to the next section 
the question of how the wide contact pads at the ends 
affect the measurement, and in the present section we 
first examine the extent to which the geometry of the 
T intersection can reduce the critical current because of 
current crowding at the sharp corners. Wc then calculate 
the boundaries of a T intersection with rounded corners 
that should prevent any reduction of the critical current. 



1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 







(a) 













-1.5-1.0-0.5 0.0 0.5 1.0 1.5 

xlW 



1.0 
0.5 
0.0 
^ -0.5 
-1.0 
-1.5 
-2.0 







(b) 













A. Sharp corners 
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The current flow in a strip of width W across the top of 
a T intersection with a voltage-contact sidebar of width 
26 is shown in Fig. 15. The conformal mapping^^ 



aw) 
aw) 



dC{w) 2b ^JW^^ 



dw 



26 



iW- —sm-\w/l3) 

TT 

2W , _i 
tanh , ^ 



(127) 



(128) 



where /3 = y^l + W^/b^, maps points in the upper half 
ti;-plane {w = u + iv) into the T-shaped region in the (- 
plane (^ = x+iy) shown in Fig. 15. The inverse mapping 
w{C,) must be obtained numerically. 

The complex potential G(^{C) is given by Eq. (7) with 
I = KjW and the inverse mapping w{Q. The corre- 
sponding complex sheet current /C^(C) = dQ(^{C,)/dQ = 
Kx{x,y) - iKy{x,y) is 



(129) 



For < y < W^, Kx{x,y) Kj as |a;| oo. The 
streamlines of the sheet current K = XK^ + yKy are 
obtained as contours of the stream function S{x,y) = 
'^Qt; {x + iy) , the imaginary part of Qi^ {() , shown in Fig. 
15(a). Note the current crowding at the inner corners of 



FIG. 15: (a) Current flow in a strip carrying total current 
KiW across the top of a T intersection, shown by the con- 
tour plot of the stream function S{x,y) — QQi^{x + iy), which 
has the values S = KjW along the top of the T {y = W) 
and S — along the underside {\x\ > b, y = 0) and the sides 
(\x\ = b, y < 0). The contours correspond to streamlines of 
the sheet-current density K, and the arrow shows the cur- 
rent direction, (b) Vortex-gcncratcd current flow, shown by 
the contour plot of the stream function Sv{xv,yv; x, y), which 
has the values Sv = along the boundaries. The contours, 
shown here for {xv,yv) = {—b + 0.1W,0.1W), correspond to 
streamlines of the vortex-generated sheet-current density Ky , 
and the arrow shows the direction of the current. The plots 
show the behavior when 6 = W/2. 



the T intersection at {x, y) = (±6, 0), where K = |/C^(C)| 
diverges. 

The critical current of the T intersection is reached 
when a vortex can be nucleated at the sharp corners in 
Fig. 15(a) or (b) at {x,y) = (±6,0). Expanding Eq. (127) 
about w = j3 (which corresponds to C, = —6) yields for 
C^^ = -b + ^e*''/^, w(C„) = p + 5w, where 



6w 



(2/3) 



1/3 ( 
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(130) 



Thus for r,^ = (—6 + (5/\/2, 5/^/2)^ where (5 ^ 6, wc can 
follow the procedure of Sec. IV to obtain the Gibbs free 
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energy, 



G = 



(t>o , /3(5\ , _ / VK\ 1/3/3(5 \ 2/3 
m ' 



27r/UoA V ^ 



(y) • ^'''^ 



Following the steps that led to Eq. (14), we obtain with 



3<5b\2/3 3</,o/3i/3 



47r2/3^oAK-/Wi/3 



5,= — ^ = 1.49^(132) 



= ^^ii, wherei?=5f^)'^' (133) 

is the reduction factor due to current crowding at one 

of the sharp inner corners, and /3 = ^/Y+JWJFp. The 
barrier height for Kj < Kc is 



(134) 



Note that the prefactor is larger than that in Eq. (15) by 
a factor of 3/2, which arises from the term proportional 
to (5^/3 in Eq. (131). 

In the limit when T4^ <C 6, /? 1, in which case Eq. 
(133) yields R « (3/2)(7r$/2VK)i/3, the same as that in 
Eq. (100) in the limit when W a and a 1. 

In the opposite limit when W ^ b, we have /3 — W/b, 
in which case Eq. (133) yields R « (3/2)«/2fe)i/3. The 
steps leading to Eq. (100) should be valid when b ^ 
but the above approximations fail for very small contact 
widths 26 ~ ^; in this case, Eq. (100) yields i? « 1. When 
2b <C ^, current-flow perturbation by the contact lead 
(the bottom of the T) is negligibly small, the current- 
crowding effect essentially disappears, and the critical 
sheet current is practically the same as in a straight long 
strip [Eq. (14)]. 



The corresponding complex sheet current JC^{C) = 
dGciO/dC, = Kx{x,y) - iKy{x,y) is 



KiWwiC) 
b^l3^-w^0' 



(136) 



For < y < W, Kx{x,y) — )• as x ^ ±00, and for 
|a;| < b, Ky{x,y) —KjW/b as y ^ —00. The stream- 
lines of the sheet current K = + yKy are obtained as 
contours of the stream function S{x,y) = ^Q(^{x + iy), 
the imaginary part of GciO^ shown in Fig. 16. Note the 
current crowding at the inner corners of the T intersec- 
tion at {x,y) = (±6,0), where K = |/C,^(<J)| diverges. 

The dashed contours in Fig. 16 for S = KjW/2 and 
3KjW/2 correspond to the optimally rounded contours 
shown in Figs. 6 and 12. As one moves along any con- 
tour under a dashed curve in Fig. 16, K = \K\ has a 
maximum near a sharp corner, but K varies monoton- 
ically when one moves along one of the dashed curves. 
Thus, for the type of current flow shown in Fig. 16, a 
patterned film in the shape of a T with rounded corners 
consisting of the region between the dashed ciirves and 
the straight line across the top of the T has the opti- 
mum shape. When ^ <^ W or b, the critical current 
will be unaffected by the bend and will be the same as 
for a long, straight film. To calculate the coordinates Xo 
and ?/o of the optimally rounded dashed curves, note that 
in the w-plane the corresponding contours are defined by 
"^oiv) = Uo + iVo = cosh ?7 -|- z sinh ?7 = sin(7r/4-|-i?7) for 
S = KiW/2 and Wo{r]) = Uo + iVo = — cosh 'q + i sinh r] = 
•^2 sin(— 7r/4 + 177) for S = 3KiW/2, where in both cases 
•q >0. Thus the coordinates of the dashed curves in the 
^-plane can be expressed with the help of Eq. (128) via 
the parametric equations 

Coiv) = Xoiv) + iVoiv) = Ciwoiv))- (137) 



B. Rounded corners 

Based upon our findings in Sees. VII B and VIII B for 
180-degree turnarounds and right-angle turns, it should 
be possible to design T intersections with rounded cor- 
ners and slightly widened strips near the intersection such 
that the critical current is not determined by vortex nu- 
cleation at the corners but rather is the same as for a 
long, straight strip. The mathematical form for a can- 
didate T intersection with rounded inner corners can be 
derived as follows. 

Let us use the same conformal mapping as in Eqs. (127) 
and (128) but consider equal current flow from the right 
and left ends of the top of the T (width W) at current 
density Kj into the bottom of the T (width 26) at cur- 
rent density KiW/b, as shown in Fig. 16. The complex 
potential describing this current flow is 



QdC) = -— Hv^HO- I]- 



(135) 



A T intersection with rounded corners as described 
above would a good choice for a patterned structure to 
avoid the problem of critical-current reduction due to 
current crowding discussed in Sec. IX A. 



X. WIDE END PAD 

To supply current to a narrow strip, it is common to 
use wide contact pads at the ends of the strip. In Sec. 
XA we show that with such a geometry the sharp cor- 
ners where the strip connects to the end pad are weak 
points where vortex nucleation caused by current crowd- 
ing occurs at a current less than the critical current of 
the long strip alone. However, in Sec. XB we calculate 
smooth boundaries of the connection between the strip 
and the end pad for which there should be no reduction 
of the critical current due to current crowding. 
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FIG. 16: Current flow in a strip carrying currents KjW in 
opposite directions merging at the top of a T intersection, 
shown by the contour plot of the stream function S{x, y) = 
"^QQix + iy) [Eq. (135)], which has the values S = 2KiW 
along the right boundary [x > b, y = and x = b, y < Q), 
S = KiW along the top of the T {y = W), and 5 = along 
the left boundary (x < —6, y = and x — —b, y < 0). 
The contours correspond to streamlines of the sheet-current 
density K, and the arrows show the current direction. The 
right dashed curve shows S = KiW/2 and the left one shows 
S = 3KiW/2. Current crowding occurs only for contours 
near the sharp corners under the dashed curves. The plot 
shows the behavior when 6 = W/2. 



A. Sharp corners 

Consider the current flow in a strip of width W (x < 0, 
—W/2 < y < W/2) connected at a; = to a wide contact 
pad of width a > W {x > 0, -a/2 < y < a/2). The 
sheet-current density is Kj for a; ^ — and KjW/a for 
x 3> a, and to describe the overall current flow we use 
the conformal mapping'^^ 



CM 

aw) 
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dw 7r(l — w'^)y/j — w ' 
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where 7 = {a^ + W'^)/{a^ — W"^), maps points in the 
upper half w-plane (w = m + iv) into the region of the 
C-plane [C, ~ x + iy) defined by < y < W/2 for a; < 
and < y < a/2 for a; > 0. The inverse mapping w(C) 
must be obtained numerically. The mapping of Eq. (139), 
shown for a ~ 2W and 7 = 5/3 in Fig. 17, corresponds to 
just the upper half of the strip and its contact pad, since 
the sheet-current density K{x, y) has mirror symmetry 
about y — 0. 

The complex potential GciO describing total current 
KjW in the x direction is given by Eq. (7) with / — 
—KjW/2 and the inverse mapping w(Q. The corre- 
sponding complex sheet current /C^(C) = dGi^{C)/dC = 



.OF 






-1.0 


-0.5 


0.0 


0.5 


1.0 






x/W 








FIG. 17: (a) Current flow in the upper half of a strip of to- 
tal width W carrying total current KiW into a wide con- 
tact strip of total width a, shown by the contour plot of 
the stream function S{x,y) = 9C/^(a; + iy), which has the 
values S = along the upper boundary and S — —KiW/2 
along the x axis. The contours correspond to streamlines of 
the sheet-current density K, and the arrow shows the cur- 
rent direction. Current crowding leading to K > Kj occurs 
for all contours above the dashed contour for S — —KiW/A. 
(b) Vortex-generated current flow, shown by the contour plot 
of the stream function Sv{xv,yv',x,y), which has the values 
Sv ~ along the boundaries. The contours, shown here 
for {xv,yv) ~ (0.05W, 0.45iy), correspond to streamlines of 
the vortex-generated sheet-current density Kv, and the ar- 
row shows the direction of the current. The plots show the 
behavior when a = 2W. 



Kx{x,y) ~ iKy{x,y) is 



(140) 



For a; < and < y < W/2, we obtain Kx{x, y) Kj as 
w{C,) — >■ —1 and X — > —00, and for x > and < y < a/2, 
we obtain Kx{x, y) Ki/2 as w(C) — > +1 and x -f 00. 
The magnitude of the sheet current, K — \JCq{C,)\, is zero 
at (a;,y) — (0, a/2), where 'w{C,) — 7, but diverges at 
(a:,y) = {Q,W/2), where w(C) 00. The streamlines 
of the sheet current K = + yKy are obtained as 
contours of the stream function S{x,y) = 55^(a; + iy), 
the imaginary part of GciOi shown in Fig. 17(a). Note 
the current crowding near the sharp corner at (a;, y) = 
(0, W/2), where K = |/C(;(C)| diverges. 

The critical current is reached when a vortex can be nu- 
cleated at (a:,y) = (0, W/2) in Fig. 17(a) or (b). To calcu- 
late it, we first need the work term W/(r„) = (/)oA/(r„) = 
(f)o[S{0,W/2) — S{xv,yv)]. However, we need this only 
very close to the corner. Expanding Eq. (138) about 
w = ±00 yields, for C = iW/2 + S/V2 - iS/V2 and 
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Thus for r„ = ((5/\/2, 14^/2 - S/V^), where (5 < VF, we 
can follow the procedure of Sec. IV and use Eq. (141) to 
obtain the Gibbs free energy, 



G = 



'"(f)-**(^)'l 



W 



1/3 / 35 N 2/3 



(7+1)' 



(142) 

Following the steps that led to Eq. (14), we obtain with 



^.,where.= 3(^y(3,. 
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is the reduction factor due to current crowding at the 
sharp corner {x,y) = (0, VK/2), and 7 + 1 = 2a? /{a? — 
W"^). The barrier height for Ki < Kc is 



47r^o A \Ki/ 



(145) 



Note that the prefactor is larger than that in Eq. (15) by 
a factor of 3/2, which arises from the term proportional 
to 52/3 in Eq. (142). 

In the limit when a/W 00, j+1 ^ 2, in which case 
Eq. (144) yields R w (3/2)(7r^/l^)i/3, the same as that 
in Eq. (100) in the limit a 1 for a strip of width W/2. 

The steps leading to Eq. (144) should be valid when 
(a—W) but the assumptions fail when {a — W) ~ ^, 
in which case (7 + 1) ^ W/S, '3> 1, and i? ~ 1. When 
(a — W) <^ ^, current-flow perturbation at the transition 
to the wider end is negligibly small, and the critical sheet 
current is practically the same as in a straight long strip 
[Eq. (14)]. 



B. Rounded corners 

We now seek the mathematical form of optimally 
rounded boundaries of a transition region from a narrow 
strip to a wide strip, such that the critical current is not 
reduced by vortex nuclcation at the corners but instead 
is the same as that of a long, straight, narrow strip. We 
already have determined in Sees. VII B and VIII B the 
shapes of optimally rounded boundaries for which there 
is no critical-current reduction due to current crowding 
around 180-degree turnarounds or right-angle turns (see 
the dashed curves in Figs. 6 and 12). Here we use a 
similar approach. 

A careful analysis of the contours for —KjW/A < S < 
in Fig. 17(a) reveals that K, the magnitude of the cur- 
rent density, has a maximum greater than Kj as one 
moves along the contour from left to right. However, for 



S = —KiW/A, the maximum occurs at a; = —00, where 
K = Ki. For contours with -KiW/2 < S < -KiW/i, 
K decreases monotonically below Kj as one moves along 
the contour from left to right. The shape of the optimal 
contour therefore can be obtained from S = —KjW/A. 
In the w;-plane this corresponds to a circular arc of ra- 
dius 1 centered at the origin, Wo{r]) = Ug + iVg = e^^ = 
COST] + isinr], where < 77 < tt. Thus the coordinates 
{xo{ri),yo{'>])) of the dashed curve in the ^-planc, the up- 
per boundary of the optimally rounded transition region, 
can be expressed with the help of Eq. (139) via the para- 
metric equations 

Uv) = Xoiv) + iVoiv) = CK(??)). (146) 

Because of the mirror symmetry about the x axis, the 
lower boundary of the optimally rounded transition re- 
gion is defined by the coordinates {xo{r]), —yoir]))- 

Note that since the rounded transition region is essen- 
tially carved out of the middle of a long strip of width W 
and an end pad of width a, far from the joint the widths 
are W/2 for the narrow strip and a/2 for the end pad. 



XI. EDGE DEFECTS 

In the above calculations, we have calculated the criti- 
cal currents of long, straight thin strips of uniform thick- 
ness or strips with straight or smooth edges interrupted 
by corners or turns. We now calculate the critical-current 
reduction when the edges of the strips are not straight or 
smooth but have defects. 



A. Semicircular notch 

Consider a long, straight strip of constant width W 
with a defect along the edge, modeled as a semicircu- 
lar notch of radius a at the edge, where a <^ W. The 
conformal mapping^** 



- 2V^ 



(2a)2 



«;(C) = C + aVC, 



{2ay 



(147) 

(148) 
(149) 



maps points in the upper half w-plane (w = u + iv) into 
the region of the C-planc (( = x + iy) above the x axis 
for |a;| > a and above the semicircle for < a. 

For a sheet-current density with the value = —Kj 
at large distances above the notch, the complex potential 



IS 



GdO = -KMC), 



(150) 



and its imaginary part is the stream function, shown in 
Fig. 18(a). Similarly, when a vortex is at the position 



23 




FIG. 18: (a) Current flow around a semicircular notch of ra- 
dius a near the edge of a strip, shown by the contour plot of 
the stream function S{x,y) = QQdx + iy), Eq. (150), which 
has the values 5 = along the boundary. The contours corre- 
spond to streamhnes of the sheet-current density K, and the 
arrow shows the current direction, (b) Vortex-generated cur- 
rent flow, shown by the contour plot of the stream function 
Sv{xv,yv', X, y), which has the values Sv ~ along the bound- 
aries. The contours, shown here for {xv,yv) = (0,1.2a), cor- 
respond to streamlines of the vortex-generated sheet-current 
density Ky, and the arrow shows the direction of the current. 
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FIG. 19: Numerical results for (a) Sc/S, and (b) R as functions 
of the ratio of the coherence length ^ to the radius a of the 
semicircular notch. Expansions in powers ol^/a are shown as 
dashed curves for 5c [Eq. (153)] and R [Eq. (154)]. Expansions 
in powers of a/^ are shown as dot-dashed curves for 5c [Eq. 
(155)] and Koc [Eq. (156)]. 



Cv = ipv, the stream function can be calculated from 
Eq. (9) with w{Q obtained from Eq. (149). The stream 
function Sy, the imaginary part of QvciCv, C)> is shown in 
Fig. 18(b). 

When a vortex is at the position = iy^ (t/„ > a) on 
the y axis, we obtain from Eqs. (10) and (11) 



G = 
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Setting i5 = y^ — a and following the steps outlined in 
Sec. IV for arbitrary values of ^/a, we obtain 5c, the 
position for which G = 0, as shown in Fig. 19(a). The 
corresponding sheet-current density is 



R, 



(152) 



where R, the reduction factor due to current crowding 
near the top of the notch, is shown in Fig. 19(b). Ana- 
lytic expansions (including only the first few terms in the 
series) for ^/o <C 1 are [dashed curves in Fig. 19] 
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Regardless of the size of ^ relative to a, vortex nucle- 
ation always occurs when the barrier height is reduced to 
zero at a distance 5c of the order of ^ from the semicir- 
cular notch. When f <SC a, the critical-current reduction 
factor is R= 1/2, which arises from current crowding at 
the top of the notch, where the sheet-current density is a 
factor of two larger than far away from the notch. How- 
ever, because of the radius-of-curvaturc effect discussed 
in Sec. VI, R approaches 1 as ^/a increases to large val- 
ues. This behavior of i? as a function of ^/a appears to 
be a general feature of edge dc^fects in thin films: When 
the linear dimensions of the edge defect are much smaller 
than ^, the suppression of the critical current is negligible 
(i?«l). 



B. Triangular notch 

Consider a long, straight strip of constant width W 
with a defect along the edge, modeled as a triangular 
notch in the shape of an isosceles triangle with base 2h 

along the edge, height a, two equal sides of length c 



and corresponding expansions for ^/a ^ 1 are [dot- Va^ + , and vertex angle Oq = 2 tan ^(6/a), where a <C 
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FIG. 20: (a) Current flow around a isosceles triangular notch 
(height a, base 26, and sides c) near the edge of a strip, 

shown by the contour plot of the stream function S{x, y) = 
^Qq{x + iy), Eq. (161), which has the values S* = along 
the boundary. The contours correspond to streamlines of 
the sheet-current density K, and the arrow shows the cur- 
rent direction, (b) Vortex-generated current flow, shown by 
the contour plot of the stream function S^{x^,yv;x, y), which 
has the values Sv = Q along the boundaries. The contours, 
shown hero for (xv,yv) = (0,0.9c), correspond to streamlines 
of the vortex-generated sheet-current density , and the ar- 
row shows the direction of the current. The plots show the 
behavior when fj, = 2/3 and the vertex angle is do = 7r/2. 



W. We define n = 1/(2 — 6o/n), which varies between 
1/2 (when 9o = 0) and 1 (when = tt). The conformal 



mapping' 
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dw " ' (1 — u;2^ 271-5 

TT /I 



Jsr(/i)=cexp - 1) 5f(;u), 



(157) 



(158) 



(159) 



C,{w) = Kjiw 2Fi(^ 



1 1 



1 1 



2^'2/z 2'2/i 



-|-zcsin 



(160) 



and 2^i(a,/3;7; z) is the hypergeometric function, maps 
points in the upper half lo-plane {w = u + iv) into the 
region of the C-plane {(^ = x + iy) above the x axis for 
I a; I > b and above the notch for \x\ < b. The inverse 
mapping w(C) must be obtained numerically. g{p) = 1 
at II = 1/2 and 1, and rises smoothly to a maximum of 
g = ^/r2(3/4) = 1.180 at M = 2/3. 

For a sheet-current density with the value Kr^ = —Kj 
at large distances above the notch, the complex potential 
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FIG. 21: The notch's critical-current reduction factor R, Eq. 
(164), versus the vertex angle 6q for ^/a = 0.01, 0.03, and 0.1, 
where a = ccos(^) = csin[Ki — 1)]. 



is 



g^iO = -KigcwiO, 



(161) 



and its imaginary part is the stream function, shown in 
Fig. 20(a). Similarly, when a vortex is at the position 
Cv = iyv, the stream function can be calculated from 
Eq. (9) with w{() obtained as the inverse of C{w), Eq. 
(160). The stream function S^, the imaginary part of 
Gv(:{Cv,C), is shown in Fig. 20(b). 

At distances along the y axis a short distance S above 
the peak at z ~ x -\riy = ia, which corresponds to w = 0, 
one can show from Eqs. (157) and (158) that w = iv, 
where v = {5/cfig)^. Following the procedure of Sec. IV, 
we obtain 



G = 



2'kijlqK 



ln(^)-0oi^,5c( 



5 

ligc) 



(162) 



Following the steps that led to Eq. (14), we obtain with 



r Sh Y_ (t>o 



5r.= 



R, where R 



H\2gcJ 



(163) 

'(164) 



is the reduction factor due to current crowding near the 
top of the notch, valid only for small values of ^/c. Figure 
21 exhibits plots of R vs 9q/ii, where 6*0 is the notch's ver- 
tex angle, and Table I exhibits the functional dependence 
of R upon ^0 and /x. The result in Eq. (164) is analogous 
to those found in Refs. 36 and 37 as the solutions of re- 
lated problems. The barrier height for Ki < is 



27r/U/UoA 



In 



(165) 



Note that the prefactor is larger than that in Eq. (15) by 
a factor of l//i, which arises from the term proportional 
to in Eq. (162). 



25 



TABLE I: The notch's critical-current reduction factor R, Eq. 
(164), as a function of the vertex angle 60 and = 1/(2 — 

eoM- 



R 




7r/4 
7r/2 
37r/4 



1/2 
4/7 
2/3 
4/5 
1 



1.414(e/a)i/2 

1.191(C/a)^''^ 
1.004(C/a)i/3 
0.876(5/0)'/^ 
1 




0.990 



1.010 



XII. THERMAL EXCITATION OVER THE 
GIBBS FREE-ENERGY BARRIER 

In the above sections we have calculated the critical 
current Ic = KcW by defining it as the current I = KjW 
at which the Gibbs free-energy barrier is reduced to zero. 
This is the critical current that would be measured in the 
limit of zero temperature, T = 0. However, because the 
Gibbs free-energy barrier height Gt is proportional to 
Hic/I), as shown in Eqs. (15), (19), (57), (101), (109), 
(134), and (145), experiments performed at finite temper- 
ature may observe switches into a state of finite voltage at 
currents less than the zcro-Gibbs-frce-cnergy-barrier crit- 
ical current Ic- Such switches can occur because nascent 
vortices may be thermally excited to the top of the free- 
energy barrier during the time of the experiment. Be- 
cause this is a statistical process, the switching currents 
can be described by a probability distribution that de- 
pends upon temperature, sample properties, and experi- 
mental conditions. 

We next derive the switching probability distribution 
applicable to experiments in which the current I through 
the superconducting device is ramped from zero to some 
maximum value at a rate r/ = dl/dt. where t is the time, 
such that / = rjt. Let Pu{t) denote the probability that 
the device is still unswitched (in the zero-voltage state) 
at time t. We assume that the probability per imit time 
of a switch into a finite-voltage state (resulting from a 
successful transit of a nascent vortex over the Gibbs free- 
energy barrier) is 



= W, I > Ic, 



{i/icr, i<ic 



(166) 
(167) 



where a; corresponds to the sample-dependent attempt 
frc^quency in Hz and is the ratio of Gb's energy scale 
to UbT: N ~ (j>Q/2TTiJ,QAkBT for a straight strip or a 
gentle curve [Eqs. (15) and (19)], N = (^oAMoAfcsT 
for a sharp 180-degree turnaround [Eq. (57)], or = 
3(/)Q/47r/ioAA:BT for sharp right-angle inner corners [Eqs. 
(109), (134), and (145)]. For the cases of interest here, 
N ^ 1. The rate of decrease of Pu{t) is given by the 
equation 



dPu/dt = -nPu. 



(168) 



FIG. 22: Switching probability distribution due to thermally 
activated barrier-climbing, P'sw/P'sw max(Ic) vs I /Ic [Eqs. 
(171) and (172)], where PL.maMc) = (iV//c) exp[-iV/(iV 
1)], shown for N = 1000 and n/ujlc. = (a) 10"^ (b) 3 x 10~^ 
(c) 10-^ (d) 3 X 10-^ (e) 10-3 {Ima. = Ic), (f) 3 X 10-^ (g) 
10-^ and (h) 3 x lO'^. 



The solution of Eqs. (166)-(168), expressed in terms of I, 
is 



= exp(^- 



UJ. 



jN+l 



exp 



{N+l)riI^)' 



Kir 



(169) 
), /> Jc.(170) 



The probability that the device has switched by the time 
the current reaches / is Pg^, = 1 — P^, and the switching 
probability distribution P'^^ = dPgw/dl = —dPu/dl is 



P' 



ri 



wl^+' 
{N+l)riI^ 

Lolc \ I Ul{I-Ic) 



OJ ( LOlc \ ( 

=;7'"H-(]vTiv7j^^H- 



ri 



I < Ic, (171) 
■),/> 4(172) 



Figure 22 shows normalized plots of P'^^ vs I/Ic for 
= 1000 and values of Ti/loIc spanning over three or- 
ders of magnitude. 

When Nrijbjlc < 1, the switching probability distri- 
bution P^^ has its maximum P'^^ at / = Imax , where 



7e(iVr,/.;/,)i/(^+i) 



(173) 



and 



PL,ma. = {N/Ima.) exp[-Ar/(7V + 1)]. (174) 



At very low temperatures N ^ 1, and for a wide 
range of values of rj/wlc we find tha,t I max 

« Ic and 

PL,ma. ~ N/elc. Since P^^dl = I, the width A/ 
of the probability distribution is approximately elc/N . 
Note that Imax = Ic when Nn/ujlc = 1- 

When Nn/bjlc > 1, the maximum of P'^^ becomes 
cusplike and remains at / = J^. For high ramp rates 
r/ ^ Lolc/N, the probability of switching at a reduced 
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current I < Ic becomes very small, and the switch- 
ing occurs chiefly for / > /c on a time scale of or- 
der 1/lu. At zero temperature, which corresponds to 
N = oo, Eqs. (171) and (172) show that no switch- 



ing occurs for I < 1^ 



j/ri at I = Ic, and 



(w/r/) exp[— (a;/r/)(/ — Ic)], with a width A7 of 



order rj/to. 



XIII. COMPARISON WITH EXPERIMENT 
USING SHARP RECTANGULAR 180-DEGREE 
TURNAROUNDS 

In 2009, Yang et al." observed that narrow hairpin 
turns in superconducting thin films exhibited reduced 
critical currents and hypothesized that the current- 
crowding effect described above might have explained 
their results, but at the time the authors were not aware 
of the work of Hagedorn and Hall, and further did not 
attempt a careful analysis.^ Yang ct al.^'* measured the 
critical currents of superconducting meander structures 
with rectangular 180-degree turnarounds like that shown 
in Fig. 11. Combining displaced mirror images of rectan- 
gular turnarounds in a two-dimensional layout (a "bous- 
trophedonic" pattern) results in a fill factor / = W/p, 
where p= W+g is the pattern period (pitch) and W and 
g are the strip and gap widths far from the turnarounds. 

The results of Sec. VII E can be used to obtain a sim- 
ple theoretical prediction for the dependence of the ratio 
of the critical current Ic{f) at fill factor / to /c(0), the 
critical current in the limit of infinite gap width (zero 
filling factor). Since in the latter case h/a 1. the de- 
sired critical current ratio obtained from 1)\ Eq. (100) 
is Ic(f)/IM = where a = ^2{h/a) - {h/af. 

However, in terms of the fill factor / we have h/a = 
(1 — /)/(! -|- /). The theory therefore predicts 



Ic{f) ^ (1 + 2/ -3/^) 

7e(0) (1 + /)V3 



2U/6 



(175) 



except for a crossover at / w 1 — ^/W into the limit 
/ ^ 1, where Sees. VII A and VII E predict 



Je(l) ^ 24/3/7r^y/6 

/c(0) i \w) ' 



(176) 



Figure 23 shows a comparison between the experimen- 
tal results and Eqs. (175) and (176). The data are the 
same as reported in Ref. 14, while the theoretical fit in- 
cludes a single free parameter, the critical current /c(0) 
in the limit of zero fill factor (/ = 0) , which was not mea- 
sured. The fitted value of /c(0) = 17 iiA was obtained 
by varying it until the fit appeared acceptable by eye. A 
few-percent variation in this fitting parameter resulted in 
a markedly imacccptablc fit. 

According to Eq. (100), since a fill factor of / = 
corresponds to the value a = 1, this critical current (17 
juA) is smaller by the factor R = (3/2)(7r^/2W)i/3 than 
the ideal critical current of a long strip of width W. With 
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FIG. 23: Comparison between experimental results for sharp 
rectangular 180-degree turnarounds^'* and the corresponding 
theoretical predictions of Eqs. (175) and (176). 



the values ^ = 7 nm and W = 90 nm^'^, R = 0.744, 
suggesting that the critical current of a long strip of width 
90 nm at the same temperature should be 23 /xA. 



XIV. DISCUSSION 

In this paper we have developed a relatively sim- 
ple, general method for systematically estimating the 
critical current of narrow, thin-film superconducting 
strip patterns with various layouts including turns and 
turnarounds. We have shown that if the latter have sharp 
inner corners, current crowding has the effect of reducing 
the critical current below that of a long, narrow super- 
conducting strip of constant width. 

Our results also have important consequences for thin- 
film superconducting single-photon detectors. ^"^^ These 
detectors carry currents close to the critical current while 
waiting for photons to arrive, but spurious dark counts 
arise when when thermally excited vortices climb over 
the Gibbs free energy barrier and cross the strip, thereby 
producing a voltage pulse. Our results suggest that the 
the frequency of dark counts in single-photon detectors is 
increased at sharp corners, where the Gibbs free energy 
barrier is reduced. 

We have made a number of simplifications to obtain 
our results, and following are various extensions and im- 
provements that could be made: 

1. Numerical solutions of the time-dependent 
Ginzburg-Landau (TDGL) equations in two dimen- 
sions, allowing for the growth of fluctuations and 
instabilities leading to vortex nucleation at the corners 
and subsequent propagation across the current-carrying 
strip, would provide more accurate values for the non- 
thermally activated critical currents. The geometries 
of interest here, for which the strip widths are much 
less than the Pearl length A — 2X^/d, provide the op- 
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portunity for considerable simplifications in the TDGL 
calculations, since the spatial variation of the order 
parameter and the current density can be calculated 
ignoring the effect of self-fields. 

2. The present calculations, which assume that the co- 
herence length ^ is much less than the strip widths, could 
be extended by relaxing this assumption. The price to be 
paid, however, is that the results for the critical current 
would become much more complicated. 

3. In calculating the self-energy, we have accounted 
only for the kinetic energy of the supercurrent circulat- 
ing around the vortex outside the vortex core. (The mag- 
netic field energy contribution is negligibly small when A 
is much larger than the strip width.) The accuracy of the 
self-energy calculation could be improved by using a vari- 
ational method to include the contributions to the vortex 
energy inside the vortex core arising from loss of conden- 
sation energy and the kinetic-energy cost of bending the 
magnitude of the order parameter."^* 

4. The calculations in this paper make use of the 
London-model assumption that the superconducting or- 
der parameter is not suppressed by the current density. 
The accuracy of the present results could be improved 
by accounting for the current-induced suppression of the 
order parameter, as done in Refs. 10 and 11. 

5. In this paper we have calculated the critical current 
as that for which the Gibbs free-energy barrier for nu- 
cleation of a vortex is reduced to zero. We believe that 
extensions of our theory to calculate the rate of thermal 
excitation of vortices over the Gibbs free energy barrier at 
slightly lower currents would confirm that the frequency 
of dark counts^"-'^^ (if time-resolved) and the thermally 
activated resistance^ (if time-averaged) are systemati- 
cally increased at sharp corners as a result of current 
crowding. 

It is important to note that in this paper we have as- 
sumed that, in stark contrast to most earlier experiments 
measuring the critical current in type-II superconductors 
with much larger transverse dimensions, bulk pinning 
plays no role whatsoever in determining the critical cur- 
rent in thin and narrow films under self-field conditions.^ 
In 2G coated conductors (thickness ~ 1 /xm and width 
^ 4 mm) with strong pinning, what normally dominates 
the critical current is bulk pinning,"^® and except in very 
rare cases^° the critical-current density due to bulk pin- 
ning is typically far below the Ginzburg-Landau depair- 
ing critical-current density.^^ However, edge pinning has 
been shown to play an increasingly dominant role in mea- 
surements using narrow strips or bridges (width < 10 
lira) to assure that the current supply is adequate to 
do the measurements. The theoretical calculations in 
the present paper are intended to apply to even thinner 
and narrower superconducting strips (thickness ~ 5 nm, 
width ^ 50-200 nm). In such films, vortices introduced 
into the strip by a large applied perpendicular magnetic 
field^^ presumably can be pinned by bulk pinning sites. 
However, if the applied field is removed, a high current 
can drive these vortices out of their pinning sites, causing 



the vortices to annihilate with their images upon exiting 
the strip. New vortices cannot enter the strip until the 
current is high enough that the vortices can surmount the 
Gibbs free-energy barrier at the other edge of the strip. 

Mirror images of the 180-degrec turnarounds discussed 
in Sec. VII can be combined in a two-dimensional lay- 
out to produce periodic meanders with the filling factor 
/ = W/p, where p = W + g is the pattern period (pitch) 
and W and g are the strip and gap widths far from the 
turnarounds. In Sec. VII B we described how to pattern a 
180-degree turnaround with an optimally rounded inner 
corner that would avoid current crowding and thereby 
maintain the critical current at the same level as that 
of an infinitely long straight strip. However, this optimal 
design yields a filling factor of only / = 1/3, because here 
W = a/2 and g = a. If it is desired to increase the filling 
factor above 1/3, our calculations indicate that this must 
come at the expense of the critical current, because, as 
shown in Sec. VII C, even with rounded inner corners the 
critical current generally will be reduced, primarily be- 
cause of unavoidable current crowding where the inner 
boundary has its minimum radius of curvature. Figure 7 
and Eq. (63) present our calculations of the best one can 
do in alleviating critical-current reduction by rounding 
the inner corner of the 180-degree turnaround. Engi- 
neering considerations rc^garding the trade-offs between 
critical currents 1^ and filling factors / will determine 
the shape of the turnaround for a specific application. 

The work we have done here in fitting this theory to 
prior work^^ suggests that the bias current in conven- 
tional SNSPDs may be limited by sharp corners. One is 
tempted to conclude from this result that appropriately 
designed devices (in which the inner corners are rounded 
to avoid current crowding) would be capable of far supe- 
rior performance. However, this presumption would be 
slightly premature-every fabrication process is slightly 
different, and the process used in the Yang paper was 
unique in its details from other processes reported in the 
literature. More careful experiments are needed in this 
case. 

Although the focus in this paper has been on supercon- 
ducting thin films, our results have relevance to the prop- 
erties of normal-metal films. The current flows shown in 
Figs. 1, 2(a), 4(a), 6(a), 8, 10, 11(a), 12(a), 13, 15(a), 
16, 17(a), 18(a), and 20(a) apply equally well to normal 
films. Current crowding at the inner corners of sharp 
bends leads to locally increased dissipation, increasing 
the electrical resistance of the strip, producing excess 
ohmic heating, and possibly increasing electromigration. 
These are all undesirable properties that could be avoided 
by choosing film patterns that optimally round the inner 
corners as discussed above. 
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of a strip. We have used this complex potential, com- 
bined with conformal mapping, to find the self-energy of 
a vortex in other strip geometries. Using the conformal 
mapping, 



((w) = ie 
w{() = ie~ 



/2a^7z/a 



(Al) 
(A2) 



we can describe the current flow surrounding the vortex 
in a thin-film sector of angular width a centered on the 
y axis (see Fig. 24) in terms of the complex potential in 
the C plane, 
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sector of opening angle a, shown by the contour plot of the 
stream function Sv{xv,yv;x,y), which has the values Sv = 
for (a;, y) along the boundaries. The contours, shown here for 
rv = {xv,yv) ~ (0, 5) and a — n/2, correspond to streamlines 
of the vortox-gcncrated sheet-current density K^, and the 
arrow shows the direction of the current. 
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i(t>o , (w{C,)-w*{C,^) 



In 



TT/ioA V w{C,) - w{Cy) 



(A3) 



whose imaginary part is the stream function 

Sv{xv,yv',x,y) = '^GvC,{Cv]C)- When the vortex is 
at (xy^yy) = (0,(5) and ^ <C (5, the circulating current 
Icirc{'''v) = 'S'u(0, (5; 0, (5 -|- ^) can be evaluated from Eq. 
(A3) by replacing the numerator in the argument of the 
logarithm by w{(v) — w*{Cy), where w{(^v) = iS'^^", and 
the denominator by (,dw{Cy)/dCv = (vrf /a)(5''/«-^ The 
resulting self-energy Eseif (ry) = 4>oIcirc{'rv)/'^ is 



Eself{5) = 



(A4) 



Appendix A: Self-energy of a Pecirl vortex in a 
sector of angle a 

Equation (8) gives the complex potential describing 
a vortex interacting with its image at the straight edge 



When the factor 2a/7r is an integer, it corresponds to the 
number of quadrants within which the vortex-generated 
current flows. Note, for example, the factor 3 in the 
self-energy term of Eq. (106) [Fig. 12(b)] and the factor 
4 in Eq. (54) [Fig. 6(b)]. 
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